Technical  Report  No.  2 


<! 


PC 

&H 


PC 


G|  W 

M 

CO 

M 
be 


o3 


00 

r> 


An  Experimental  Study 
of  Wind  Velocity  Profiles 
over  a  Wavy  Surface 

by 

Francis  J.  Merceret 


ONR  N000I4-69-A0407 


iTgr^nMEKT  A 

~T  A  ;  -v  rivblic  release; 

Ap'«  —a  '  •  . 


O  D  C 

__  ,, — N  r — i  pp  '  1  j 

n\E£P'  "V 

Ijj  HIM  5  V312 


Dii-.nlDuaon 


Unlimited 


Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

Springfield.  Vs  221  SI 


jun  » 
■gHSETie 


COLLEGE  OF  MARINE  STUDIES 

UNIVERSITY  OF  DELAWARE 

Newark,  Delaware  .  / 


AN  EXPERIMENTAL  STUDY 


OF  WIND  VELOCITY  PROFILES 
OVER  A  WAVY  SURFACE 


by 


Francis  J.  Merceret 


Lecturer,  College  of  Marine  Studies 
University  of  Delaware 
Newark,  Delaware  19711 


Details  of  iHcstnvionr;  in 
this  docenKir'l  ir.sy  be  batter 
siiniicd  on  initroliehe 


The  text  of  this  document  is  also  submitted  to  the  Department  of 
Earth  and  Planetary  Sciences,  The  Johns  Hopkins  University, 
Baltimore,  Maryland,  as  a  dissertation  (1972,  unpublished)  for  the 
degree  of  Doctor  of  Philosophy. 


Supported  by  ONR  Contract  Number  N00014-69-A-0407 
Research  in  the  Coastal  and  Oceanic  Environment 


Publication  2MS055  of  the 
College  of  Marine  Studies 
at  the  University  of  Delaware 


Newark,  Delaware  19711 
April  1972 


Unclassified _ 

_  Security  Ch>^siHc_otion _ _ _ ____________ _ 

P*^5  DOCUMENT  CONTROL  DATA  RiD 

.  SVcur/fy  t  f«>.iiUraHi>ri  <»/  tl.tlo,  body  of  nbstrnct  and  /ndtM/ritf  onnoMtlon  niunt  be  entered  when  the  ovcrmtl  repot!  In  clam  tilled) 


1  O  M  l  O  l  N  A  TINO  AC  TIVITV  (i'orporate  tuilhot)  2».  REPORT  SECURITY  CLASSIFICATION 

|  Unclassified 

University  of  Delaware 


]  «EPORT  TITLE 

An  Experimental  Study  of  Wind  Velocity  Profiles  Over  a  Wavy  Surface 


4  OESCfliP  TIVE  NOTE*  (Type  ut  report  find  inc/iisfvr  dm  ten) 

Technical  heport 


S  AuTHOPiSi  (Fire!  middle  inltiel,  Imet  name) 

Francis  J.  Merceret 


A  REPO*  T  O  A  TT 

April,  1972 


7*.  TOTAL  NO.  OP  PACES  7b.  NO-  OF  PEPS 

113  29 


If.  CONTRACT  OR  GRANT  NO- 

NOOOU|-69-A-OU07 

b.  PRO  JCC  T  NO. 

9».  ORIGINATOR’?1  REPORT  NUMBE  R  (S) 

Technical  Report  No.  2 

c. 

d. 

9b.  other  report  noISI  (Any  other  number *  thmt  may  be  meelgned 
thin  report) 

CMS  No.  2MS06$ 

10.  DISTRIBUTION  STATEMENT 

Distribution  of  this  document  is  ■ 

unlimited. 

11-  SUPPLEMENTARY  NOTES 

12.  SPONSORING  Ml  LI  TAR  Y  ACTIVITY 

Geography  Branch 

Office  of  Naval  Research 

Arlington,  Virginia  22217 

■  13.  ABSTRACT 

Wind  velocity  profiles  and  turbulence  measurements  have  been  obtained 
over  a  plow  corrugated  field  using  hot-film  anu  cup  anemometers  at 
sampling  heights  covering  from  3  cm  to  1*  m.  The  profiles,  based  on 
averaging  times  from  $  to  1 $  minutes,  fit  the  logarithmic  model  with 
correlation  coefficients  greater  than  0.9  and  with  many  greater  than  0.97. 
The  data  support  the  conjecture  of  Ruggles  (1970)  that  large  and 
sporadic  fluctuations  in  the  roughness  length  and  friction  velocity 
observed  over  water  are  due  to  the  generation  of  wind  waves  on  the  surface. 
Additional  information  regarding  the  response,  of  the  profile  parameters 
to  rapidly  changing  conditions  is  inferred  from  studies  of  profiles 
inside  and  outside  of  intermittent  bursts  of  turbulent  activity.  The 
results  indicate  that  mechanisms  which  involve  steady  state  logarith¬ 
mic  mean  velocity  profiles  in  the  wind  above  a  water  surface  are  not 
suitable  to  explain  the  generation  of  waves,  since  the  data  show  the 
profile  parameters  to  be  markedly  affected  by  the  generation  process. 


,1473  < 


PAGE  I  ) 


Unclassified 


S/N  0)0  1  ■  807-  680  1 


Security  Classification 


Waves 


Turbulence 


Wind  profiles 
Air-sea  interaction 


DD  Vov.,1473  <back> 

(PAGE  2) 


Unclassified 

Security  Classification 


ABSTRACT 


Wind  velocity  profiles  and  turbulence  measurements  have  been 
obtained  over  a  plow  corrugated  field  using  hot- film  and  cup 
anemometers  at  sampling  heights  covering  from  3  cm  to  U  m. 

The  profiles,  based  on  averaging  times  from  5  to  15  minutes, 
fit  the  logarithmic  model  with  correlation  coefficients 
greater  than  0.9  and  with  many  greater  than  0.97*  The  data 
support  the  conjecture  of  Ruggles  (1970)  that  large  and  spora¬ 
dic  fluctuations  in  the  roughness  length  and  friction  velocity 
observed  over  water  are  due  to  the  generation  of  wind  waves  on 
the  surface.  Additional  information  regarding  the  response  of 
the  profile  parameters  to  rapidly  changing  conditions  is  in¬ 
ferred  from  studies  of  profiles  inside  and  outside  of  inter¬ 
mittent  bursts  of  turbulent  activity.  The  results  indicate 
that  mechanisms  which  involve  steady  state  logarithmic  mean 
velocity  profiles  in  the  wind  above  a  water  surface  are  not 
suitable  to  explain  the  generation  of  waves,  since  the  data 
show  the  profile  parameters  to  be  markedly  affected  by  the 


generation  process. 
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I .  INTRODUCTION 


The  experiment  described  in  this  paper  was  conducted  over  land, 

t 

yet  it  is  directed  towards  understanding  the  generation  of  sur¬ 
face  water  waves  by  the  turbulent  wind.  Because  it  is  not 
immediately  obvious  that  one  ought  investigate  wave  generation 

in  a  place  where  no  waves  are  being  generated,  it  seems  appro- 

\ 

priate  to  explain  why  and  how  the  decision  to  make  this 
particular  investigation  evolved.'. 

In  a  series  of  papers,  Miles  (1957 »  1959,  I960,  1962)' 

\  l 

presented  a  model  which  has  become  the  basis  for  mlich  of  cur¬ 
rent  research  into  the  generation  of  surface  waves  by  the  wind. 
Expositions  on  and  interpretations  of  the  Miles  model  may  be 
i’ouhd  in  Lighthill  (1962),  Kinsman  (1965),  and  Phillips  (1966) 
among  others,  while  a  summary  of  the  model  and  a  discussion  of 
its  place  in)  over-all  wind  wave  theory  is  given  in  Shemdin 
and  Hsu  (1966).  because  others  have  written  so  extensively 
on  the  subject,  I  will  describe  here  only  those  features  of 
the  theory  pertaining  to  the  purposes  and  results  of  this  in¬ 
vestigation. 

Miles  calculated  the  energy  and  momentum  transfer  from 
a  laminar  shear  flow  of  prescribed  velocity  profile  to  a  train 
of  water  waves  after  the  wind  field  had  been  modified  by  the 
presence  of  the  waves.  He  assumed  the  velocity  profile  in  the 
air  was  steady  and  that ' the  mean  and  periodic  components  in  the 
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wind  were  separately  two  dimensional  though  not  necessarily 
parallel.  The  waves  were  assumed  not  to  modify  the  mean  wind 

field.  Turbulence  was  included  only  insofar  as  the  mean  pro- 

\ 

files  suggested  for  use  in  the  model  were  those  common  to 
turbulent  rather  than  laminar  flows.  The  logarithmic  profile 
was  considered  the  most  applicable  to  prediction  of  ocean  wave 
development. 

Laboratory  and  field  investigations  showed  that  Miles' 
theory,  while  better  than  its  competitors  t  consistently  under¬ 
estimated  the  growth  rates  of  waves.  Shemdin  (1969)  compiled 
several  observers'  results  and  compared  them  with  his  own.  He 
concluded,  in  agreement  with  Phillips  (1966)  and  others,  that 
the  omission  of  turbulent  fluctuations  and  second  order  pro¬ 
ducts  of  the  wave  perturbations  was  a  major  reason  for  the 
failure  of  the  theory  to  estimate  correctly  the  rate  of  energy 
transfer. 

Experimental  work  continued,  while  other  investigators 
turned  to  the  task  of  including  the  neglected  effects  in  a 
more  general  theory.  The  introduction  of  turbulence  into  a 
model  always  brings  with  it  the  problem  of  closure.  The 
averaged  equations  for  any  quantity  in  a  turbulent  flow  al¬ 
ways  contain  correlations  of  one  order  higher,  i.e.,  the, 
equations  for  the  mean  flow  contain  the  double  velocity  and 
pressure  velocity  correlations,  those  for  the  double  velocity 
correlations  contain  triple  velocity  correlations  and  sc 
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ad  infinitem.  Most  of  the  expansions  on  Miles'  vork  thus  settled 
down  to  introducing  turbulence  into  the  equations ,  choosing  a 
closure  scheme,  and  solving  the  resulting  equulions  to  see  if 
anything  believable  resulted. 

In  order  to  solve  a  set  of  differential  equations,  one 
must  have  boundary  conditions.  For  typical  models  such  as 
those  derived  by  Hussain  and  Beynolds  (1970),  Manton  (l97l), 
and  Merceret  (unpublished)  it  is  necessary  to  specify  boundary 
conditions  at  a  surface  and  at  infinity  or  another  surface  for 
both  the  mean  flow  and  the  perturbations.  For  solid  walls, 
such  as  those  used  by  Hussain  in  his  experimental  study,  this 
poses  no  problem.  For  application  to  water  waves  under  a  semi¬ 
infinite  atmosphere,  however,  it  is  impossible  to  establish 
sensible  boundary  conditions  without  knowing  the  answer  before 
one  begins  the  calculation.  The  models  are  created  in  order  to 
determine  the  transfer  of  energy  and  momentum  to  the  waves.  If 
one  knows  the  stresses  at  the  surface  exactly,  the  transfer  can 
be  computed.  If  not,  then  the  boundary  conditions  cannot  be 
specified. 

A  more  profitable  approach  seemed  to  me  to  be  to  seek 
relationships  among  the  variables  in  the  region  near  the  surface 
which  would  enable  me  to  reduce  the  number  of  free  variables 
and  use  these  and  the  kinematic  boundary  condition  to  provide 
sufficient  conditions  at  the  surface.  Like  the  closure  schemes 
used  in  turbulence,  this  would  require  a  combination  of 
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observation,  dimensional  analysis,  and  physical  reasoning,  and, 
like  them,  would  give  an  approximate  answer. 

I  had  decided  to  make  a  systematic  and  careful  field 
study  of  the  lowest  meter  of  the  atmospheric  boundary  layer 
using  hot-film  anemometry  to  see  whether  there  existed  a 
detectable  region  of  transition  from  a  horizontally  homogenous 
turbulent  flow  with  a  logarithmic  profile  to  a  region  of  dif¬ 
ferent  characteristics,  possibly  a  laminar  sublayer.  1  would 
also  investigate  the  details  of  the  transition  region  if  I 
found  it.  If  I  could  find  a  laminar  sublayer  and  relate  the 
flow  in  it  to  that  outside  I  could  use  the  laminar  boundary 
conditions  of  Miles  and  our  new  information  to  solve  the  ex¬ 
tended  Miles  model.  Even  if  no  laminar  sublayer  existed,  I 
might  find  clues  to  useful  relations  for  what  did  exist. 

The  study  would  begin  on  land  and  go  to  sea  only  if 
something  useful  resulted  from  the  land-based  work.  Time  and 
money  were  limited  and  we  could  not  afford  the  risk  and  expense 
involved  in  working  at  sea  until  I  was  sure  I  had  a  chance  of 
generating  knowledge  worth  the  cost  of  its  acquisition. 

That  was  what  I  planned  to  do.  What  was  actually  done 
was  something  more  important. 

While  I  was  waiting  for  my  equipment  to  be  delivered, 
we  received  a  visit  from  Professor  Erik  Mollo-Christensen  of 
MIT,  who  presented  a  lecture  at  a  Civil  Engineering  seminar. 

His  lecture  cast  doubt  on  the  whole  framework  of  assumptions 
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underlying  the  Miles  theory  and  in  private  discussion  he  sug¬ 
gested  the  experiments  I  finally  decided  to  perform. 

Much  of  what  Dr.  Mollo-Christensen  presented  of  the  work 
of  his  group  has  not  yet  been  published.  The  material  here  is 
from  their  published  work  together  with  what  Dr.  Mollo-Christensen 
suggested  to  me. 

If  one  goes  to  sea  in  search  of  logarithmic  velocity  pro- 
f iles  to  substitute  into  Miles '  equations ,  one  finds  them. 
Disturbingly,  though,  the  roughness  parameter  zQ  Beems  to  depend 
strongly  on  the  wind  speed  and  on  the  properties  of  the  surface 
in  a  way  which  Ruggles  (1970)  suggests  is  a  function  of  the 
generation  mechanism.  The  data  of  Barger  et  al.  (1970)  support 
Ruggles'  hypothesis,  as  do  those  of  Mollo-Christensen  (1970). 

A  more  detailed  discussion  of  theae  data  will  be  presented  in 
Chapter  II. 

If  the  parameters  of  the  profile  depend  on  the  generation 
of  waves  at  the  surface,  one  cannot  assume  a  profile  prior  to 
computing  the  growth  rate  of  the  waves  as  required  for  Miles' 
theory  and  those  models  currently  derived  from  it.  The  search 
for  proper  boundary  conditions  must  wait  until  we've  determined 
whether  we  yet  have  a  proper  set  of  averaged  field  equations. 

If  Ruggles  is  right,  the  best  available  theories  explaining  the 
generation  of  waves  by  turbulent  wind  are  wrong  in  principle 
and  effort  should  be  directed  from  expanding  on  them  to  devising 
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a  more  correct,  fully  interacting  model. 

A  major  reason  that  the  significance  of  the  MIT  work  has 
been  less  appreciated  than  one  might  expect  may  be  that  there 
has  been  nothing  with  which  to  compare  it.  The  variation  of 
z0  with  wind  speed  might  be  a  function  of  the  geometry  of  the 
boundary  rather  than  its  mobility.  Wind  tunnel  measurements 
are  not  adequate  because  of  the  large  difference  in  Reynolds 
number  between  laboratory  and  field  situations.  The  bursts  of 
high  frequency  turbulence  reported  by  Mollo-Christensen  (1970) 
associated  with  rapid  changes  in  profile  parameters  might  be 
causing  the  changes  rather  than  occurrences  at  the  surface. 

When  I  discussed  my  proposed  program  with  Dr.  Mollo- 
Christensen,  he  agreed  that  our  instrumentation  and  site  were 
excellent  to  provide  the  supporting  data  to  confirm  or  confute 
his  group's  work.  Dr.  Mollo-Christensen  agreed  that  this 
would  be  a  more  fruitful  use  of  our  limited  resources  so  I 
altered  my  plans  accordingly. 
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II.  THE  STUDY 

The  possibility  that  the  mean  wind  profile  is  not  independent 
of  the  wave  field  is  not  original  with  Ruggles ,  but  earlier 
papers  imply  a  general  belief  that  any  variation  of  roughness 
length  will  be  a  smooth  one  as  the  wind  speed  increases, 
Priestley  (1959) i  Kraus  ( 196T )  >  or  at  least  piecewise  smooth 
with  a  single  region  of  transition  where  boundary  layer  sepa¬ 
ration  over  the  waves  begins,  Bole  and  Hsu  (1967).  The  data 
of  Ruggles  (1970)  indicate  the  reality  to  be  more  complicated. 

Ruggles  fitted  data  from  nearly  300  runs  taken  at  sea 
with  fetches  greater  than  26  miles  to 

lnz  =  —rU(z)  +  lnz 
u*  o 


where 

z  is  the  height  above  the  reference  level, 
k  is  Von  K&rm£n's  constant, 

U  is  the  mean  wind  speed, 
u*  is  the  friction  velocity,  and 
z0  is  the  roughness  length. 

Table  1  displays  the  surprisingly  large  number  of  close  fits 
he  obtained.  For  these  data  the  roughness  length  and  friction 
velocity  behaved  in  an  unexpected  manner.  Fig.  1,  Fig.  2. 

Note  that  the  large  ranges  of  u*  for  a  given  wind  speed  at 


Correlation  coefficient  Number  of  samples 


1.00 

86 

0.980-1.00 

106 

0.950-0.980 

59 

0.920-0.950 

26 

0.890-0.920 

11 

0.860-0.890 

2 

0.830-0.860 

3 

0.800-0.830 

2 

0.7UO-0.800 

0 

0.700-0. 7^0 

h 

Table  1.  Correlation  with  logarithmic  profile 
from  Ruggles  (1970) 


.Lj 


from  Ruggles  (1970)  from  Ruggles  (1970) 


10  meters,  U1(j,  occur  where  the  peaks  in  the  zQ  plot  occur. 

The  multiple  peaks  disappear  in  the  presence  of  an  oil  slick, 
Fig.  3,  Barger  et  at.  (1970)-  The  large  ranges  indicated  for 
u*  do  not  represent  values  of  u*  distributed  over  the  whole 
range.  Rather,  the  values  of  u“  cluster  at  the  extremes  of 
the  range,  and  the  value  appearing  in  a  particular  profile 
depends  on  the  ratio  of  presence  to  absence  of  high  frequency 
bursts  in  the  turbulence,  Mollo-Christensen  (1971). 

The  MIT  group  interpret  these  data  as  showing  wind- 

speed  dependent  momentum  interchange  between  the  wind  field 

and  the  wave  field.  When  large  transfer  occurs  the  surface 

of  the  sea  roughens  as  capillary  and  small  gravity  waves  are 

generated  on  the  backs  of  the  waves  constituting  the  swell. 

Thus  increases.  At  the  same  time  the  transfer  of  momentum 
o 

changes  the  values  of  the  Reynolds  stresses,  thus  changing 
u*.  The  burst  structure  for  the  high  frequency  turbulence 
coincides  with  a  similar  phenomenon  on  the  water  surface  where 
small  wavelets  are  seen  to  be  generated  in  patches  which  sud¬ 
denly  appear  and  disappear,  Kinsman  (personal  communication). 

The  burst  structure  has  been  observed  in  the  atmosphere 
over  land,  Haugen,  Kaimal,  and  Bradley  (1971),  and  in  the 
laboratory,  Kim,  Kline,  and  Reynolds  (1971)  but  the  variation 
of  roughness  length  with  wind  speed  has  been  observed  only  over 
yielding  surfaces  such  as  the  sea  or  very  long  pliant  grass, 
Priestley  (195 9)>  The  roughness  length  over  flat  land  seems 
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independent  of  wind  speed,  Priestley  (1959),  Lumley  and 
Panofsky  (1964).  The  value  of  zq  as  a  function  of  wind  speed 
has  not  been  systematically  investigated  over  a  periodic  fixed 
surface,  as  far  as  we  ore  aware,  outside  of  the  laboratory; 
thus  it  is  still  possible  that  Ruggles'  data  may  be  indicative 
of  a  shear  instability  in  the  wind  flow  rather  than  of  occur¬ 
rences  at  the  surface. 

The  present  experiment  is  designed  to  test  the  hypothesis 
that  the  surface  configuration  is  alone  responsible  for  the 
value  of  zq,  even  over  a  corrugated  surface,  and  to  examine  the 
burst  structure  as  it  relates  to  the  value  of  u*  in  logarithmic 
profiles. 

In  order  to  get  the  strongest  possible  connection  with 
Ruggles'  work,  an  experimenter  should  use  similar  instruments 
in  a  similar  configuration  and  process  his  data  in  the  same 
manner  as  Ruggles.  Of  course,  additional  information  beyond 
what  Ruggles  obtained  should  be  used  when  it  is  available  to 
illuminate  the  problem  and  care  should  be  taken  to  ensure  that 
the  results  are  not  a  fiction  of  the  instruments  or  data 
handling  methods.  The  MIT  experiments  used  light-interrupting 
cup  anemometers  for  mean  wind  speeds  at  higher  levels  and 
hot-film  anemometers  at  lower  levels  for  mean  winds  and  tur¬ 
bulence  levels.  We  had  ordered  the  same  kinds  of  instruments, 
which  will  be  described  in  a  later  section.  The  field  site 
should  be  one  periodically  corrugated  with  a  long  unobstructed 
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fetch.  Such  a  site  was  available  to  us  and  will  be  described 
in  the  next  section. 

In  the  experiment  we  would  have  to  do  the  following: 

1.  Obtain  velocity  profiles 

2.  Select  those  which  are  logarithmic  with  a 
correlation  of  0.9  or  better 

3.  For  these  profiles  plot  zQ  va.  Uio  and  note 
the  presence  or  absence  of  peaks  such  as 
Ruggles  found 

4.  Plot  u*  va.  U,Q  and  look  for  regions  of  high 
variance  and,  if  they  exist,  correlate  them 
with  such  peaks  as  show  up  in  the  zQ  plot 

5.  Look  for  bursts  of  high  frequency  activity  at  . 
each  level  and  determine  whether  they  occur 
simultaneously,  sequentially,  randomly,  or  not 
at  all 

6.  If  bursts  occur  nearly  simultaneously,  reduce 
portions  of  the  record  during  bursts  and  the 
remainder,  separately  fitting  them  to  logarith¬ 
mic  profiles. 

The  goodness  of  fit  should  yield  information  on  the  rate  of 
adjustment  of  the  profile  while  the  values  for  zq  and  u*  may 
be  used  to  determine  whether  u*  correlates  with  the  burst 

structure  and  whether  z  is  a  function  of  the  burst  mechanism. 

o 

If  the  bursts  occur  sequentially  with  the  same  ratio  of  time 
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occupied  by  bursts  to  that  of  the  total  record  obtaining  at 
each  level  then  this  method  would  still  be  used,  but  additional 
information  would  be  gained  from  studying  the  lag  time  between 
levels.  If  the  bursts  occur  randomly  at  each  level  or  not  at 
all,  it  would  tend  to  indicate  the  phenomenon  is  not  of  suf¬ 
ficient  depth,  or  frequency  respectively,  to  be  a  major 
determining  factor  in  the  shape  of  the  profile. 

If  the  burst  structure  is  observed,  and  if  z  is  not  a 
function  of  the  wind  speed  or  the  burst  structure  but  u*  is 
dependent  on  both,  then  the  speculations  of  the  MIT  group  are 
well  founded  and  we  need  to  reconsider  our  acceptance  of  the 
prescribed-profile  models  accordingly.  If  otherwise,  then 
some  investigation  into  other  mechanisms  for  the  variation  of 
the  profile  parameters  would  be  indicated. 
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III .  '  THE  SITE  '  \  ! 

l 

\ 

I  preferred  a  field  site  in  the  Delavare  Valley  as  close^  to 

I  : 

Newark  as  could  be  found  because  my  laboratory  and  data  analysis 

facilities  were  in  Newark,  Delaware  at  The  University  of 

\  '  \ 

Delaware's  College  of  Marine  Studies.  \  I  needed  a  periodically 

corrugated,  relatively  unvegetated  surface  with  a  long  unob¬ 
structed  fetch  in  the  direction  of  the  predominant  wind.  The  j 

■  \ 

experiment  would  require  a  wind  relatively  steady  in  direction 

and  mean  magnitude  over  periods  of  the  order  of  tens  of  minutes. 

Thus  the  time  of  the  year  for  the  work  and  the  alignment  of 

i 

the  fetch  were  determined  by  the  climatology  of  the iDelaware 

i 

Valley  region. 

Th^  winds  steadiest  in  direction  and  of  sufficient 

I 

speed  to  cover  the  range  of  interest  occur  in  Delaware  from 

\ 

late  October  to  early  April  with  the  months  February  and  March  ; 
providing  the  best  conditions,  l .  .ther  (196Q).  During  these 
months  and  during  early  April,  the  prevailing  wind  is  from 
the  west  tc  the  north  with  northwest  winds  occurring  nearly 
30 %  of  the  time.  Wind  roses  for  March  and  April  for  the 
station  closest  to  our  site  are  shown  in  Fig.  4,  Because 
our  equipment  arrived  in  late  February  we  operated  only  during 
the  period  mid-March  to  mid-April.  The  actual  winds  were  quite 
satisfactory  during  this  time,  being  from  the  desired  direction 
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over  a  suitable  range  of  speeds. 

The  vicinity  of  The  University  of  Delaware  abounds  in 
small  airports  and  plowed  fields.  A  plowed  field  is  suitably 
corrugated  while  an  airport  usually  has  an  adequate  fetch  in 
the  direction  of  the  prevailing  wind.  We  were  fortunate  to 
find  a  plowed  field  adjacent  to  a  well  kept  but  little  used 
grass  airstrip,  the  owners  of  each  of  which  were  willing  to 
tolerate  a  clutter  of  equipment,  a  tower,  and  interloping 
scientists  on  their  property.  The  airport's  windsock  was  a 
bonus.  From  it  we  could  estimate  the  steadiness  of  the  wind 
direction.  Figure  5  shows  the  site  which  is  located  near  the 
Delaware-Maryland  border  in  Elkton,  Maryland.  On  Fig.  6, 
note  the  location  of  the  site  and  of  Delaware  City  for  which 
the  wind  roses  have  been  presented. 

In  the  direction  of  the  predominant  wind,  northwest, 
the  fetch  is  unobstructed  and  the  surface  almost  level  for 
more  than  1  km.  There  were  no  obstructions  in  the  semicircle 
from  SW  through  NW  to  NE  closer  than  about  900  meters  except 
for  several  small  trees  and  a  number  of  two  story  farm  buildings 
due  west  at  about  550  meters,  and  a  stand  of  trees  slightly 
east  of  north  at  about  480  meters.  Downwind  obstructions  in¬ 
cluded  a  two  story  house  and  some  one  story  aircraft  hangers 
about  360  meters  away  to  the  SSE  and  SE,  respectively,  and  a 
stand  of  trees  due  east  at  about  500  meters.  Otherwise,  the 
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downwind  fetch  was  unobstructed  for  at  least  700  meters  in  all 
directions. 

The  site  is  plowed  from  the  farm  buildings  to  the  edge 
of  the  airstrip.  The  plow  contours  have  the  approximate  form 
of  an  unsymmetrical  square  wave  with  an  amplitude  of  10  cm  and 
a  wavelength  of  90  cm,  Fig.  7.  The  furrow  crests  run  nearly 
east-west.  Corn3talks  ranging  from  a  centimeter  to  perhaps 
30  cm  high  and  1-3  cm  in  diameter  punctuate  the  surface  every 
meter  or  so,  but  many  of  them  have  been  flattened  against  the 
ground  so  that  those  remaining  form  a  random  and  relatively 
sparse  pattern.  The  soil  was  firm  even  after  heavy  rain,  and 
water  did  not  appear  to  collect  anywhere  in  the  field  during 
the  two  months  we  used  the  area;  nor  was  the  field  at  any  time 
snow  or  ice  covered. 

The  site  is  ideal  for  this  research.  Its  proximity  to 
our  laboratory  made  it  possible  to  test,  repair,  and  reinstall 
equipment  during  a  day's  run.  It  enabled  us  to  operate  on 
very  short  notice  when  the  weather  was  right.  The  steady  wind 
and  long  unobstructed  fetch  aligned  nicely,  while  the  nature 
and  contours  of  the  surface  were  as  if  made  to  order. 


>ntour  of  a  typical  furrow 
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IV.  THE  INSTRUMENT  SYSTEM  AND  ITS  FIELD  INSTALLATION 

The  wind  instruments  were  cup  anemometers  and  hot-film  anemometers. 
The  cup  anemometers,  including  data  registers,  were  C.W.  Thornthwaite 
Associates  Wind  Profile  Recording  System  Model  106.  The  anemometer 
cups  are  of  light-weight  plastic  reinforced  with  aluminum  rings 
while  stainless  steel  tubing  forms  the  spokes  and  the  hub  is  of 
aluminum.  The  entire  cup  assembly  weighs  only  7  grams.  The 
anemometer  shaft  is  made  of  hardened  stainless  steel.  Miniature 
ball  bearings  are  vised.  A  shutter  mounted  on  the  shaft  inter¬ 
rupts  a  light  beam  from  a  3-volt  grain-of-wheat  lamp  housed  in 
the  base  of  the  anemometer  tee  and  breaks  a  beam  of  light 
focused  on  a  photocell  once  for  every  revolution  of  the  cup 
assembly.  After  amplification  the  output  of  thiB  photocell 
operates  an  electromechanical  counter  directly.  No  relays  are 
used  in  the  Wind  Profile  Recording  Systems.  One  count  registers 
on  the  counter  for  each  revolution  of  the  cup  assembly.  The 
counters  are  mounted  on  the  front  panel  of  the  cabinet  for  ease 
in  reading.  The  anemometers  in  our  set  were  factory  calibrated 
and  matched  to  within  +0.2$,  The  calibration  curve  for  the 
system  is  shown  in  Fig.  8.  The  specifications  supplied  by  the 
manufacturer  and  confirmed  at  Delaware  are  shown  in  Table  2. 

The  response  of  the  cup  anemometers  to  tilt  is  essentially 
flat  within  1+5°  of  horizontal,  see  Fig.  9,  Thornthwaite,  et  al. 
(1959),  and  the  system  acts  like  a  low  pass  filter  for  small 


WIND  SPEED  RANGE: 


ilk 

0-32  mph,  0  -  1U50  cm/sec 
STARTING  SPEED: 

leas  than  0.20  mph  (8.9^  cm/aee) 

DISTANCE  CONSTANT: 

83.3  cm  in  horizontal  flow 

READOUT  RESOLUTION: 

0.1  cm/sec  (l  horn'  profile) 

POWER  SUPPLY: 

1  -  12  v  storage  battery 
1  -  1.3  v  mercury  cell 

INDICATION: 

electromechanical  register 
ANEMOMETER  CUPS: 

conical,  plastic,  2  in  diameter;  assembly  weight, 

7  grams  (0.25  oz) 

TRANSMITTER  HOUSING: 

chrome  plated  brass  mounted  on  3/8  in  o.d.  arms  12  in 
long;  total  weight  of  transmitter  156  grams  (5.5  oz); 
upright  cup  assembly  support,  3/8  in  o.d.;  height 
Ik  cm;  maximum  interference  dimension  viewed  by  the 
wind,  1/2  in 

CABINET: 

9"K  x  8-1/ 8"W  x  12"D,  3.63  kg  (8  lbs) 

Table  2-  Manufacturer's  specifications  for  cup  anemometer  system 
supplied  by  C.W.  Thornthwaite  Associates 


IOC 
80 
60 
40 
cr.  20 
0 

-20 
-40 

-60  -30  0  30  60  90  -60  -  30  0  30  60  90 

Angle  from  horizontal 


Fig.  9.  Directional  response  of  cup  anemometers 
from  Thornthwaite  et  at.  (1959) 
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scale  fluctuations,  see  Fig.  10,  Thornthwaite ,  Superior,  and 
Mather  (1961}.  For  our  system,  the  3  db  response  point  occur¬ 
red  at  about  4  meters  with  essentially  full  response  for  scales 
larger  than  20  meters. 

The  hot-film  sensors  were  DISA  55A81  quartz-coated  wedge 
shaped  probes.  DISA  55E05  battery  operated  constant  temperature 
anemometers  and  55D15  portable  linearizers  were  used.  Films 
were  preferred  to  wires  because  they  are  more  rugged  and  relatively 
insensitive  to  weather,  Merceret  (1970).  We  chose  to  linearize 
the  system  for  two  reasons.  The  range  of  winds  expected  and  the 
turbulence  level  anticipated  were  large  enough  that  lineariza¬ 
tion  should  significantly  contribute  to  accuracy,  and  the  ability 
to  suppress  the  zero-wind  bridge  voltage  would  allow  better  use 
of  the  limited  dynamic  range  of  the  recording  system. 

The  probes  were  operated  at  about  200  C.  At  this  tempera¬ 
ture,  in  wind  fields  of  the  kind  expected,  the  frequency  response 
of  the  system  is  essentially  flat  to  beyond  200  Hz,  with  the 
3  db  point  well  beyond  1  kHz,  Merceret  (1969),  more  than  ade¬ 
quate  for  the  proposed  measurements,  see  Fig.  11.  The  directional 
response  of  the  films  is  not  so  forgiving  as  that  for  cups, 

Merceret  (1968),  and  a  steady  wind  and  careful  alignment  of  the 
probes  were  necessary  for  high  accuracy.  Factory  specifications 
for  the  anemometers  and  linearizers  are  presented  in  Tables  3 
and  4.  The  specifications  for  the  anemometers  were  confirmed, 
Merceret  (1968),  (1969),  while  those  for  the  linearizers  were 
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TYPES  OF  OPERATION: 


PROBES : 


10:1  bridge  ratio  and  1:1  bridge  ratio 


•  hct-wire  or  hot- film  types 
note:  hot-film  probes  normally  require  an  external 
resistor 


PROBE  RESISTANCE  RANGE: 

built-in  bridge  arm:  2.U  to  13.0  ohms  in  10  steps 
external  bridge  arm:  1  to  50  ohms 

NOISE  LEVEL: 


approximately  0.8  mVrms  s.t  zero  flow  velocity  and 

\  i 

1  50  kHz  bandwidth  (the  noise  .level  is  proportional 

to  the  bandwidth) 

OUTPUT  VOLTAGE: 

approximately  1  to  7  volts  at  full  flow  velocity 
OUTPUT  IMPEDANCE: 

approximately  500  ohms 


FREQUENCY  RANGE: 

10:1  bridge  ratio:  0-50  kHz.  The  upper  frequency 
limit  depends  on  the  conditions  of  measurement  and 
on  rhe  type  of  probe  used 

1:1  bridge  ratio:  0  -  100  kHz.  The  upper  frequency 
limit  depends  on  the  accuracy  of  the  compensating 
cable,  on  the  conditions  of  measurement,  and  on  the 
type  of  probe  used 

POWER  SUPPLY: 

built-in,  two  HELLESENS  Type  K-10  9-volt  batteries 
or  equivalent  and  six  Type  BA-30  1.5-volt  dry  cells 
or  equivalent,  for  instance  HELLESEKS  Type  VII-36 

Table  3.  Manufacturer's  specifications  for  hot-film  anemometer  system 

supplied  by  DISA  S  &  B  (continued  on  page  30) 


30 

I 

i  I 

i 

Table  3.  continued  from  page  2 9  ^ 

EXTERNAL  POWER  SUPPLY i 

extra  bridge  pover  may  be  'supplied  (maximum  20 
volts  DC) 

AMBIENT  TEMPERATURE: 

0  to  +1*5  C  • 

CIRCUIT: 

l 

fully  transistorized  (silicon  transistors) 

DIMENSIONS: 

62  mm  high,  212  mm  wide,  220  mm  .deep 
WEIGHT:  ) 


2.1*  kg  inclusive  of  batteries 
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BASIC  TRANSFER  FUNCTION: 

m  m 

V  =  K(V  -  V  )  and  V  =  K  V 
out  in  ino  out  in 

EXPONENT  (ra): 

2  -  3  -  It  -  5  -  6  (±5%) 

INPUT  RESISTANCE  (R  =  V  /I  ): 

in  in  in 

850S  to  5  kN 

MAXIMUM  INPUT  VOLTAGE: 

30  volts 

OUTPUT  VOLTAGE  RANGE: 

maximum  output  voltage  better  than  3*5  volts 
DC  output  voltage  can  be  measured  with  built-in  meter 
meter  range:  0-3  volts 
meter  accuracy:  2$  of  full  scale 

OUTPUT  IMPEDANCE: 

approximately  10 

maximum  load:  1  kfi  (or  maximum  output  current:  it  mA) 
OUTPUT  VOLTAGE  CHANGE  PER  DEGREE  CENTIGRADE  (AKB.  TEMP.): 


Typical  Maximum 


0.2$ 

per 

C 

0.5$ 

per 

C 

(•at. 

m  = 

2  and  3) 

0.3$ 

per 

C 

0.6$ 

per 

c 

(at 

m  = 

10 

0.U$ 

per 

C 

0.8$ 

per 

c 

(at 

m  = 

5) 

0.5$ 

per 

C 

1.0$ 

per 

c 

(at 

m  = 

6) 

OUTPUT  VOLTAGE  STABILITY: 

AVqu^/10  hours  at  constant  temperature 

maximum  0.05  volt  with  battery  supply- 

maximum  0.001  volt  with  external  stabilized  power 

supply 

Table  4.  Manufacturer ’ s  specifications  for  linearizers 
supplied  by  DISA  S  &  B  (continued  on  page  32) 
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Table  1* .  continued  from  page  31 

TYPICAL  OUTPUT  NOISE: 

1  mV  at  full  bandwidth 
rms 

0.2  mV  at  10  kHz 
rms 

DIMENSIONS  (HMD)  IN  MILLIMETERS: 

6 2  x  212  x  280 

WEIGHT: 

1.8  kg 

POWER  SUPPLY: 

built-in:  twelve  HELLESENS  VII  28  1.5-volt  dry  cells 
or  equivalent  (e.g.  MATSUSHITA  ELECTRIC  UM3D,  Mallory 
Radio  Penlight  M15R) 

AMBIENT  TEMPERATURE  RANGE: 

0  to  +1*5  C 

TYPICAL  UPPER  FREQUENCY  LIMIT  (3  db):  m  =  3: 


33 


confirmed  during  this  study  except  as  will  be  noted  in  Chapter  V. 

The  wind  instruments  were  mounted  several  ways.  The  cup 
anemometers  were  always  attached  to  an  8-meter  tri-sectioned  mast 
of  two  inch  aluminum  irrigation  pipe  guyed  so  as  to  avoid  inter¬ 
ference  from  the  wakes  of  the  guy  wires.  The  attachment  was 
made  using  crossarms  of  3 A  in  square  aluminum  stock  and  U  bolts, 
as  is  conventional  in  mounting  TV  antennas.  Each  cup  was  about 
eighteen  inches  from  the  mast  and  positioned  to  the  side  of  the 
upwind-downwind  line  to  avoid  aerodynamic  interference  from 
the  mast.  Fig.  12. 

We  had  four  cup  anemometers  in  the  array,  the  top  three 
of  which  were  at  H  m,  2  m,  and  1  m  above  the  top  of  a  furrow 
except  during  the  instrument  checkout  run  during  which  the 
four  cups  were  at  7*9  a,  3.9  a,  1.9  a,  and  90  cm.  The  bottom 
cup  was  placed  at  the  same  height  as  the  top  film  on  every 
run  to  cross-check  the  instruments.  This  height  was  either 
50  cm  or  25  cm  depending  on  the  run. 

The  hot-film  units  were  on  separate  mounts  described  in 
the  next  paragraph.  There  were  three  working  channels  available, 
the  top  one  of  which  was  always  at  50  cm  or  25  cm  above  the 
crest  of  a  furrow.  The  remaining  two  probes  were  then  located 
at  two  heights  selected  from  the  following:  with  a  top  of 
50  cm;  25  cm,  12.5  cm,  and  6  cm  or  with  a  top  of  25  cm;  12.5  cm, 

6  cm,  3  cm,  and  1  cm.  None  of  the  data  taken  at  the  1  cm 


Fig.  12.  Cup  anemometer  installation 
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position  were  recoverable  but  all  of  the  data  taken  at  greater 
heights  were  usable. 

The  hot-film  anemometers  were  mounted  on  18  in  steel 
shafts  attached  to  aluminum  blocks  which  fitted  onto  1/2  in 
wooden  dowelB.  The  dowels  provided  electrical  isolation  of 
each  probe.  These  dowels  were  attached  to  the  tower  cross- 
arms  for  instrument  cross  checks  and  thereafter  mounted  in 
the  ground  be3ide  the  large  tower  as  probe  supports.  Fig.  13. 

In  order  to  have  a  measure  of  the  temperature  distri¬ 
bution  in  the  air  near  the  surface,  we  installed  two  mercury- 
in  glass  meteorological  thermometers  with  radiation  shields 
at  heights  from  0.20  to  2  m  above  the  surface.  These  were 
read  to  the  nearest  0.01  C  at  the  beginning  and  end  of  each 
run. 

The  thermometer  readings  and  cup  anemometer  counts  were 
recorded  on  prepared  forms  with  the  time  and  date  of  the  begin¬ 
ning  and  end  of  each  run.  The  data  from  the  hot-film  systems 
were  recorded  on  a  Lockheed  model  417  analog  tape  recorder  and 
the  date  and  time  marked  on  each  reel. 

We  used  the  Model  417  in  a  4-channel  FM  record-playback 
configuration  using  a  Lockheed  storage  cell  for  power.  The 
tape  speed  most  convenient  for  our  work  was  7-1/2  ips  which 
gave  us  flat  response  from  D.C.  to  well  above  1  kHz  while 
allowing  enough  tape  for  30  minutes  of  data  plus  another  5 
minutes  of  calibration  signals  on  each  reel.  The  manufacturer's 
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specifications  are  presented  in  Table  5.  A  three  to  one  voltage 
divider  of  precision  resistors  was  used  to  step  down  the  voltage 
from  the  linear! zers  to  a  value  compatible  with  the  limits  of 
the  recorder.  This  enabled  us  to  run  the  55D15*8  at  full  gain 
so  there  would  be  less  probability  of  misadjusting  the  gain 
control  than  if  we  used  an  intermediate  setting,  and  it  allowed 
us  to  utilize  the  full  dynamic  range  of  both  instruments. 
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1.  GENERAL 

TAPE  SPEEDS:  7-1/2,  15,  and  30  ips 

TAPE  DIMENSIONS:  up  to  3150  feet  of  1/2  in  polyester-backed 

recording  tape  can  be  accommodated  on  a  7  in 
reel 

RUNNING  TIME:  45  minutes  at  7-1/2  ips  for  l800  feet 

START  AND  STOP  TIME:  acceleration  to  7-1/2  ips,  5  sec  maximum 

FLUTTER:  the  cumulative  peak-to-peak  flutter  as  measured  95# 
of  the  time,  will  not  exceed  0.8#  at  7-1/2  ips 


2.  FM  SYSTEM 

RMS  3IGNAL-TO-RMS  NOISE:  40  db  measured  across  passband  of 

recorder  at  7-1/2,  15,  and  30  ips 

HARMONIC  DISTORTION:  1.5#  total  harmonic  distortion  measured 

at  1/2  cut-off  frequency 

A. C./D.C.  LINEARITY:  1#  of  full  scale  at  all  speeds 

DRIFT:  1.0#  of  full  scale  from  1/2  minute  to  1  hour,  at  constant 
temperature  of  68  F.  0.3#  Of  full  scale  from  5  minutes 
to  1  hour,  at  constant  temperature  of  68  F 

TEMPERATURE  STABILITY:  0.05#  per  degree  F  from  32  to  120  F 

INPUT  SENSITIVITY:  ±1.4  volt  (1.0  volt  rms )  for  ±40#  deviation 

INPUT  IMPEDANCE:  20,000  ohms  minimum  over  the  passband 

OUTPUT  VOLTAGE:  ±1.4  peak,  volts  (1.0  volt  rms)  for  full 
carrier  deviation 

OUTPUT  IMPEDANCE:  1000  ohms  maximum  over  the  passband 
OUTPUT  LOAD:  10  K  ohms  minimum 


Table  5-  Manufacturer's  specifications  for  tape  recorder 
supplied  by  Lockheed  (continued  on  page  39) 
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Table  5*  continued  from  page  38 
MAXIMUM  DEVIATION:  ±40* 

v 

CENTER  FREQUENCY  AT  7-1/2  IPS:  13,500  Hz 
FREQUENCY  RESPONSE  AT  7-1/2  IPS:  0  to  2,500  Hz 

3-  ENVIRONMENTAL  AND  PHYSICAL  SPECIFICATIONS 

OPERATING  TEMPERATURE  RANGE:  0  to  1*5  C 

HUMIDITY:  up  to  99%  (without  condensation) 

ORIENTATION:  operates  in  any  position 

SIZE,  STANDARD  CASE:  13-15/16"  x  15-3/16"  *  6-3/8",  exclusive 

of  handle 

WEIGHT,  STANDARD  CASE:  29  pounds,  including  b elf-contained 

battery 

POWER  REQUIREMENT:  approximately  13  watts  at  a  nominal 

17.5  volts 


V.  THE  CALIBRATION  AND  USE  OF  THE  SYSTEM 


We  checked  the  factory  calibration  of  the  cup  anemometers  and 
calibrated  the  hot-film  systems  in  the  12  x  18  inch  test  section 
of  a  lov  turbulence  wind  tunnel  at  The  Johns  Hopkins  University, 
see  Kellogg  (1965),  using  a  pitot-static  tube  and  liquid  micro- 
manometer  as  the  reference  instrument  for  velocity. 

The  manufacturer's  specifications  for  the  cup  anemometers 
were  correct  to  within  about  1%  of  the  indicated  flow  at  velo¬ 
cities  above  the  starting  velocity  as  presented  in  Fig.  8. 

The  hot-film  systems  were  tested  in  the  configuration 
to  be  used  in  the  field  with  12  v  tractor  batteries  supplying 
power  to  the  probes,  and  initially  at  the  field  overheat 
temperatures.  Later  tests  were  made  at  a  variety  of  over¬ 
heat  temperatures.  With  the  no- flow  output  voltage  suppressed 
and  the  linearizer  controls  adjusted  as  they  would  be  in  the 
field,  repeatable  curves  such  as  those  shown  in  Fig.  lh  were 
obtained  for  each  system.  Deviations  from  linearity  were  at 
most  3/S  of  the  indicated  flow  at  any  point  on  the  curve. 

These  curves  seemed  to  be  independent  of  the  operating  tem¬ 
perature  of  the  probe  for  any  given  system,  a  desirable 
feature  but  one  which  required  explanation. 

The  explanation  was  discovered  during  an  investigation 
of  another  aspect  of  the  calibration  of  the  system.  If  the 
zero  suppression  was  not  done  precisely,  the  slope  and 


Fig.  1^4.  Hot-film  calibration  curve  for  zero  suppressed  system 


Fig.  15.  Hot-film  calibration  curve  for  unzeroed  system 
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Intercept  of  the  curve  changed,  Fig.  15.  The  change  in  inter¬ 
cept  was  to  be  expected.  The  change  in  slope  was  an  unpleasant 
surprise  which  indicated  that  the  gain  of  the  system  was  a 
function  of  the  zero  suppress  control  position. 

The  instruction  manual  for  the  55D15  linear! zer  gave  no 
indication  that  the  gain  would  be  a  function  of  the  zero  sup¬ 
pression  control  position.  Therefore,  I  decided  to  check  the 
units  against  the  diagram  in  the  manual.  They  agreed.  I  then 
derived  the  transfer  characteristics  of  the  unit  from  the 
circuit  diagram.  The  result  explained  what  I  had  observed. 

The  actual  transfer  function  for  the  55D15  is 

Eout  *  [C  (KEt»  -  v]""her» 

K  is  a  constant  determined  by  the  zero  suppress  control, 

VQ  is  a  constant  determined  by  the  min/max$  control, 

C  is  a  constant  determined  by  the  gain  control,  and 

m  is  a  constant  determined  by  the  exponent  control. 

The  derivations  of  this  expression  and  the  conclusions  following 
are  contained  in  a  College  of  Marine  Studies  Technical  Report, 
The  University  of  Delaware,  presently  in  press.  Because  copy¬ 
righted  material  is  involved  it  will  not  be  reproduced  here. 

If  one  substitutes  King's  Law  in  the  form 


into  the  transfer 

\ 

I 

Eout  =  «*<>>“  [<*  +  f  "  x]  m  for  a  properly 

zeroed  system.  Thus  the  gain  for  a  zeroed  system  is  independent 
of  both  the  overheat  temperature  AT  and  of  the  position  of  the 
zero-suppress  control  so  long  as  the,  system  is  balanced. 

For  a  given  overheat,  the  gain  is  clearly  a  function  of 
K  under  unbalanced  conditions.  But  the  function  is  now  known. 
Thus',  a  correction  can  be  applied  to  date  taken  while  the  control 
is  misset  if  the  resulting  offset  is  not  too  large,  Merceret 
(in  press).  The'  correction  consist^  of  a  shift  and  a  scaling, 
thus  \ 

E  t  -  4  (E  .  -  E„)  where  \ 

corrected  G  out  o' 

Eq  is  the  shift  voltage  equal  to  20U^.(U  =  0)  and  G  is  the  gain 

of  the  system  relative  to  the  gain  of  a  properly  suppressed 

system,  which  function  may  be  determined  by  computation  or  by 

measurement  for  each  unit  as  a  function  of  the  value  of  EQ. 

Because  we  had  recorded  E  in  the  field  we  were  able  to  correct 

o 

for  the  error  in  those  runs  for  which  E  was  not  exactly  zero. 

o 

Details  of  the  method  are  presented  in  the  afore  mentioned 


E  =  1(1  f-,UN)  A  At] 
anem  L  A  J 


function  for  the  lineariser,  one  gets 


lilt 

\ 

l 

CMS  Technical  Report.  Of  'the  eleven  reels  of  data  we  collected, 
six  required  no  correction,  three  were  corrected,  and  two  were 
discarded  as  being  beyond  the  range  for  which  reliable  cor’.’ections 
could  be  made. 

The  final  step  in  the  laboratory  checkout  of  the  instru¬ 
mentation  was  to  calibrate  the  film  systems  ana  cup  anemometers 
against  each  other  over  the  speed  range  of  1  -  10  m/s.  The 
systems  agreed  within  ±2%  of  the  indicated  velocity. 

We  also  tested  the  tape  recorder-  in  the  laboratory 

I 

After  an  initial  problem  was  solved  which  required  returning 

.  \ 

the  unit  to  the  factory  fo~  repair,  the  machine  performed 

within  specifications  Under  test.  Nonetheless,  in  case  the 

recorder  was  adversely  affected  by  the  temperature  or  voltage 

fluctuations  likely  in  the  field,  we  took  the  precaution  of 

\ 

placing  record  system  calibration  marks  on  the  tape.  A  switch- 
box  containing  the  voltage  divider  permitted  each  channel  to  be 
selectively  connected  to  a  3*1  v  reference  source  or  effectively 
grounded1  by  the  divider  network  before  signals  were  applied, 
see  Fig.  l6.  If  the  carrier  frequency  or  gain  of  the  recording 
system  varied,  it  wouM  show  up  as  a  zero  o:  reference  voltage 
shift  respectively  at  playback,  and  a  suitable  correction 
could  be  applied  to  the  lata.  In  no  case  on  the  field  tapes 
did  such  a  shift  actually  occur. 


VI.  THE  RECOVERING  AND  THE  PROCESSING  OF  THE  DATA 


The  thenaometric  data  and  those  from  the  cup  anemometers  were 
in  numerical  form  as  taken  and  could  he  plotted  or  fitted  to 
profiles  numerically  without  processing  except  for  converting 
cup  counts  to  velocities  by  a  simple  multiplication  in  each 
case. 

The  data  from  the  hot-film  systems  were  recorded  on 
tape  as  analog  signals  convertible  to  voltage  linearly  propor¬ 
tional  to  sensed  wind  speed  as  a  continuous  function  of  time. 
The  constant  of  proportionality  for  each  system  was  determined 
from  the  calibration  curves.  It  was  necessary,  therefore,  to 
process  the  electrical  signal  to  determine  the  mean  voltage 
for  each  run  to  obtain  profile  values ,  and  to  look  at  the 
energy  contained  in  the  signal  above  20  Hz  to  examine  the 
burst  structure.  While  we  had  originally  intended  to  use  elec¬ 
tronic  analog  processing  to  obtain  the  information  we  required 
and  to  produce  turbulent  one-dimensional  energy  spectra  as 
well  to  determine  as  much  as  possible  about  the  structure  of 
the  flow,  we  were  unable  to  carry  out  these  plans  for  financial 
reasons.  Moreover  we  found  that  the  rented  recorder  had  to  be 
returned  within  two  weeks  of  the  completion  of  the  field  work. 
Thus  we  would  lose  our  playback  facilities. 

To  avoid  the  unpleasantness  of  having  reels  of  valuable 
but  inaccessible  data  we  devised  another  method  of  processing 


the  information  which  would  give  us  all  that  we  needed  and 
most  of  what  we  wanted  except  spectra.  We  placed  a  rush  order 
for  two  Keathkit  Model  EU20B  chart  recorders  with  variable 
speed  driveB,  manufacturer's  specifications  for  which  are 
presented  in  Tables  6  and  ?•  The  data  would  be  transferred 
from  magnetic  tape  to  charts  using  the  ftistest  available  chart 
speed  of  5  sec/inch,  the  low  frequency  signal  components  on 
one  chart  and  the  high  frequency  components  on  the  other.  The 
separation  used  a  low  pass  filter  and  a  high  pass  filter  with 
their  inputs  in  parallel,  their  outputs  being  routed  to  the 
low  frequency  information  storing  recorder  and  a  rectifier- 
filter  system  respectively.  The  output  of  the  rectifier-filter 
system  went  to  the  high  frequency  information  storing  recorder. 
A  circuit  diagram  of  the  system  is  shown  in  Fig.  17*  The 
high  pass  filter  was  a  General  Radio  Model  1952  Universal 
Filter  operated  in  the  high  pass  mode  with  a  cutoff  frequency 
setting  of  20  Hz.  Manufacturer's  specifications  and  rejection 
curves  are  shown  in  Table  8  and  Fig.  16,  respectively.  The 
low  pass  filters  between  the  tape  and  the  low  frequency  chart 
recorder,  and  between  the  rectifiers  and  the  high  frequency 
chart  recorder  were  designed  to  severely  attenuate  signals 
above  10  Hz  to  avoid  causing  any  instabilities  in  the  chart 
recorders.  The  chart  recorders  have  60  Hz  choppers  in  the 
servo-system  and  the  system  can  begin  to  "hunt"  wildly  if 
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CHART  PAPER:  grid  width:  10  inches;  length:  120  foot  roll; 
markings:  0  -  100 

CHART  SPAN:  five  fixed  ranges:  10,  25,  50,  100,  and  250  mV; 

plus  a  vernier  adjustment  between  each  fixed 
range,  which  provides  adjustment  of  span  to  any 
value  between  250  and  3.3  mV.  There  is  also  an 
external  position  for  special  plug-in  fixed  ranges 

BALANCING  TIME:  0.1  second  per  inch,  1  second  full  scale 
(10  inches) 

INPUT  CIRCUIT:  self-balancing  potentiometric  on  all  ranges 

INPUT  RESISTANCE:  essentially  infinite  at  null.  Approximately 
500  kfl  off  null 

LINE  FREQUENCY  REJECTION:  80  db  in  common  mode 

FLOATING  INPUT:  100  v  DC  maximum  above  ground.  70  VAC  rms 
maximum  across  the  negative  input  terminal 
and  GND  terminal 

OVERALL  ERROR:  less  than  1 #  of  full  scale  from  10  to  250  mV 

DEAD  ZONE:  less  than  0.5#  of  full  scale  from  10  to  250  mV 
LINEARITY:  less  than  0.5#  full  scale 
MAXIMUM  RECOMMENDED  SOURCE  RESISTANCE:  50  k?2 
REFERENCE  SOURCE:  zener  regulated  supply 
POWER  REQUIREMENTS:  120  volts,  60  Hz;  45  watts 
DIMENSIONS:  l4"W  x  8-3/4"H  x  13-1/4"D 


Table  6.  Manufacturer's  specifications  for  chart  recorder 


supplied  by  Heath  Company 
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CHART  SPEEDS: 

5  seconds  to  2  hours  per  inch  in  21  speeds  when 
using  the  internal  chart  drive  signal.  Any  chart 
speed  up  to  5  seconds  per  inch  when  using  an 
external  chart  drive  signal 

CHART  ADVANCEMENT  STEPS: 

600  steps  per  inch;  2  steps  per  cycle  of  the  signal 
applied  to  the  motor  drive  circuit.  The  motor  will 
accurately  step  from  0  to  1+0  cps  and  will  run 
accurately  at  60  cps 

TRANSISTOR  COMPLEMENT: 

27  -  2N3393:  Schmitt  triggers  binaries,  and 
emitter  follower.  2  -  2N31+16:  motor  drive 

DIODE  COMPLEMENT: 

29  -  1N191  crystal  diode:  twenty-four  binary  and 
five  feedback  diodes.  1  -  1N751  zener  diode:  D.C. 
level  adjuster.  2  -  750  mA  100  v  PIV  silicon  diodes  : 
supply  rectifiers 


CONTROLS : 


21-position  chart  speed  switch.  External- internal 
signal  switch 

POWER  REQUIREMENTS: 

105-125  volts  AC,  60  cps,  30  watts 


Table  7- 


Manufacturer's  specifications  for  chart  drive 
supplied  by  Heath  Company 
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FREQUENCY  RANGE 

CUT-OFF  FREQUENCIES.-  adjustable  h  Hz  to  60  kHz  in  four  ranges 

PASS-BAND  LIMITS:  low-frequency  response  to  D.C.  (approximately 
0.7  Hz  with  A.C.  input  coupling)  in  LOW  PASS 
and  BAND  REJECT  modes .  High-frequency  response 
uniform  ±0.5  db  to  300  kHz  in  HIGH  PASS  and  BAND 
REJECT  modes 

CONTROLS:  log  frequency-dial  calibration;  accuracy  ±2%  of  cut¬ 
off  frequency  (at  3-db  points) 


FILTERS 

FILTER  CHARACTERISTICS:  filters  are  fourth-order  (four-pole) 

Chebyshev  approximations  to  ideal 
magnitude  response.  The  nominal  pass- 
band  ripple  is  ±0.1  db  (±0.2  db  max); 
nominal  attenuation  at  the  calibrated 
cut-off  frequency  is  3  db;  initial 
attenuation  rate  is  30  db  per  octave. 
Attenuation  at  twice  or  at  one-half 
the  selected  frequency  as  applicable, 
is  at  least  30  db 

TUNING  MODES:  switch  selected,  LOW  PASS,  HIGH  PASS,  BAND  PASS, 
and  BAND  REJECT 


INPUT 


GAIN:  0  or  -20  db,  switch  selected.  Accuracy  of  gain  is  ±1  db, 
of  20  db  attenuator  is  ±0.2  db 

IMPEDANCE :  100  kfi 

COUPLING:  A.C.  or  D.C..  Lower  cut-off  frequency  (3  db  down) 
for  A.C.  coupling  is  about  0.7  Hz 


Table  8.  Manufacturer's  specifications  for  bandpass  filter 
supplied  by  General  Radio  (continued  on  page  52) 
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Table  8.  continued  from  page  51 


GENERAL 

OUTPUT:  600  (2  impedance 

NOISE:  100  nV  in  an  effective  bandwidth  of  50  kHz 
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Fig.  18.  Frequeni 


60  Hz  energy  reaches  the  amplifiers.  This  filtering  also 
smoothed  the  record  somewhat  and  made  it  easier  to  read  with¬ 
out  removing  much  of  the  detail.  The  chart  pens  would  each 
respond  substantially  to  changes  with  time  scales  longer  than 
200  ms.  At  time  scales  of  800  ms  or  larger  essentially  com¬ 
plete  response  was  obtained. 

A  provision  to  place  simultaneous  marks  on  the  high  and 
low  frequency  charts  allowed  for  the  possibility  of  synchro¬ 
nizing  the  charts  even  if  the  chart  speeds  differed  slightly. 
The  precaution  proved  to  he  unnecessary,  the  recorders  always 
agreeing  within  2  seconds  after  a  30  min  run,  and  within  one 
second  from  any  intermediate  mark  to  the  closer  end  of  the 
chart  where  signal  marks  denoted  the  beginning  and  end  of 
each  tape.  Since  the  minimum  practically  resolvable  time 
difference  between  events  on  the  charts  was  about  1  second 
this  meant  the  chart  synchronization  was  for  practical  pur¬ 
poses  perfect. 

We  could  now  use  graphical  averaging  techniques  on  our 
paper  chart  record.  We  had  replaced  the  expense  of  electronic 
analysis  with  the  labor  of  manual  date  reduction. 

The  pairs  of  charts — one  high  frequency  information 
chart  and  one  low-frequency  information  chart  per  channel  from 
the  tape — were  filed  according  to  data  and  run  number  and 
marked  with  information  about  the  system  and  height  from  which 


55 


its  record  came.  Because  we  would  sample  the  mean  velocity 
selectively  if  bursts  were  occurring  in  the  high  frequency 
energy,  the  high  frequency  charts  of  each  run  were  analyzed 
first.  It  is  necessary  here  to  preview  the  results  in  order 
to  explain  subsequent  processing  of  the  data.  Bursts  were 
observed  at  all  levels  simultaneously  in  each  run  in  nearly 
every  case  where  a  burst  was  seen  in  any  record.  Figure  19 
shows  a  typical  burst,  while  Fig.  20  shows  bursts  occurring 
simultaneously  at  each  height  with  the  charts  time  synchro¬ 
nized  using  the  calibration  and  end  of  record  marks. 

Each  set  of  high  frequency  charts  containing  records 
from  three  heights  during  a  given  run  was  time  synchronized 
using  the  spikes  at  the  beginning  and  end  of  recording,  see 
Fig.  21,  and  examined  for  bursts.  Judgements  were  made  as 
to  where  bursts  appeared  in  each  set  of  records  and  the  region 
of  burst  activity  marked  on  each  chart,  see  Fig.  22.  These 
marks  were  transferred  to  the  low  frequency  records  after 
they  had  been  time-mark  synchronized  with  the  high  frequency 
graphs.  Thus  the  mean  velocity  record  could  be  selectively 
sampled. 

The  bursts  were  sufficiently  marked  that  when  I 
asked  visitors  to  my  office  to  pick  them  out,  they  repeatedly 
identified  them  and  marked  the  leading  and  trailing  edges 
within  a  second  or  two  of  my  judgement;  thus  the  method  was 
a  consistent  and  simple  one  to  use.  Any  algorithmic  system 
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such  as  marking  a  burst  whenever  the  trace  exceeded  a  certain 
level  involves  decisions  as  to  where  to  place  the  marker  level 
and  thus  the  human  judgement  is  not  removed,  it’s  just  better 
concealed.  Moreover,  there  were  trends  in  the  data  which 
would  cause  difficulty  in  consistently  applying  any  algorithm 
to  select  bursts.  Since  we  were  looking  for  a  structure  of 
a  particular  and  readily  distinguishable  kind  superimposed 
on  trends  of  other  kinds,  the  trained  human  intellect  seemed 
the  tool  most  likely  to  produce  meaningful  results . 

Having  selected  which  parts  of  each  record  were  to  be 
averaged  together  it  was  necessary  to  choose  a  procedure  for 
averaging  the  data  from  the  charts.  The  first  method  we 
considered  was  the  use  of  a  planimeter  to  integrate  the  trace 
but  repeated  trials  showed  variations  of  as  much  as  5#  in  the 
value  obtained.  We  decided  to  see  what  happened  if  we  just 
"eye-balled"  a  straight  line  through  short  sections  of  the 
graph  and  computed  average  values  from  the  line  segments. 

The  results  were  pleasantly  surprising:  the  answers  were 
always  within  ±2%  of  the  mean  of  the  planimeter  readings  and 
the  technique  was  much  quicker.  The  average  for  each  seg¬ 
ment  was  taken  to  be  the  value  at  center  of  the  straight  line 
defining  the  segment  and  the  length  oi  the  segment  was  used 
as  a  weighting  factor  in  computing  the  overall  average  on  a 
Wang  Model  500  calculator.  The  accuracy  of  the  method  is 
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confirmed  by  the  results  of  the  instrument  checkout  run,  which 
are  reported  in  Chapter  VII. 

In  addition  to  the  mean  velocity  information  we  desired 
a  measure  of  the  energy  above  20  Hz  within  and  without  the 
bursts.  Because  we  could  not  calibrate  the  detector  network 
and  filter  system  for  the  unknown  spectral  shape  of  the  in¬ 
coming  signal,  lacking  the  necessary  equipment,  the  measure 
would  have  to  be  a  relative  one  and  would  be  that  only  if  the 
spectra  above  20  Hz  maintained  approximately  the  same  shape 
regardless  of  amplitude.  In  the  inertial  subrange  this  would 
be  the  case. 

If  the  utility  of  the  relative  measure  be  granted, 
obtaining  a  mean  square  still  is  somewhat  more  difficult  than 
obtaining  a  mean.  As  with  the  means,  I  examined  mechanical 
methods  of  integrating  the  squares  of  the  values  of  the  curves 
and  tested  various  graphical  techniques.  I  came  to  the  con¬ 
clusion  that  short  horizontal  line  segments,  much  shorter 
than  those  used  for  mean  values,  should  be  employed  and  their 
position  squared  and  averaged  with  the  appropriate  weighting. 
If  the  value  of  the  line  in  any  segment  was  £,  and  the  graph's 
fluctuations  from  that  value  within  the  segment  5'  ,  then 

.2  2  2 
5  =  £  +  V 

where  is  the  true  mean  square  value  of  the  graph  in  the 
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segment.  Since  £  is  a  constant,  £  =  £  and  thus,  If  the  seg- 

~~~2  — 2  TT  —2 

ment  is  made  small  enough  that  £'  <<  £  ,  then  £  ~  £  .  The 

segments  were  picked  accordingly. 

The  Wang  Model  500,  a  programmable  desk  computer-cal¬ 
culator,  was  used  for  the  arithmetical  analysis  of  the  data. 
When  geometrical  averages  were  wanted,  a  program  for  arith¬ 
metical  averages  was  used  entering  the  logarithm  of  each  datum. 
The  antilogarithm  of  this  mean  is  the  geometric  mean.  Applying 
"exponentiation"  to  the  standard  deviation  of  the  mean  of  the 
logarithms  yields  a  value  which  pertains  to  the  deviation  in 
ratio  rather  than  the  deviation  in  value  from  the  mean. 

Because  the  directional  responses  of  cup  anemometers 
and  hot-film  anemometers  are  quite  different,  a  directionally 
varying  wind  could  cause  the  two  kinds  of  instruments  to  give 
different  results  for  the  mean  wind  at  a  given  height  and 
thus  "warp"  the  profile.  Our  method  of  dealing  with  this 
possibility,  for  the  lack  of  anything  better,  is  as  follows. 

If  we  assume  the  wind  vectors  are  parallel  at  all  heights 
sensed  by  our  instruments  then  the  ratio  of 


U( measured  cup) 
U(measured  film) 


will  be  the  same  for  any  height. 


It  will  be  greater  than  one  since  the  hot-films  respond  less 
than  fully  to  winds  at  an  angle  to  their  direction  of  maximum 
sensitivity  (for  which  the  calibrations  apply)  while  the 


cups  are  omni-directional  in  the  horizontal.  Since  we  had  a 
cup  and  a  film  together  at  one  height  we  computed  the  ratio 
and  multiplied  it  into  the  velocities  obtained  from  the  lower 
films  and  based  the  results  on  cup-equivalent  velocities. 

The  distribution  of  values  obtained  is  displayed  in  Table  9* 
Their  relative  reliability  was  considered  in  reviewing  the 
data  but  the  profiles  did  not  seem  to  be  poorer  when  larger 
factors  were  involved. 
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VII.  THE  RESULTS 


One  of  my  Immediate  concerns  following  the  completion  of  the 
field  work  was  for  the  accuracy  and  precision  which  I  could 
expect  after  the  whole  process  of  reducing  the  data  had  "been 
completed.  The  errors  of  analysis  usually  compound  the  errors 
of  measurement  and  recording.  How  good  would  the  final  result 
be? 

Before  we  began  profile  measurements  we  made  three 
checkout  runs  with  the  hot-films  mounted  side  by  side  with 
a  cup  anemometer  as  a  reference  at  90  cm  during  winds  typical 
of  those  encountered  during  the  profile  measurements.  The 
results  indicate  we  may  have  considerable  confidence  in  the 
overall  system.  In  Table  10  are  displayed  the  ratios  of  the 
mean  velocity  values  obtained  from  the  graphic  averaging  pro¬ 
cedure  for  hot- film  systems  A,  B,  and  C  to  that  obtained  from 
the  reference  cup  anemometer  in  each  run.  Of  the  nine  entries, 
five  show  an  overall  error  smaller  than  or  equal  to  3%,  and 
only  one  exceeds  7%.  These  data  were  taken  during  inter¬ 
mittent  snow  showers  and  sunlight  which  presented  the  worst 
operating  environment  in  which  we  had  to  work  and  with  the 
instruments  over  different  parts  of  the  furrows,  thus  inducing 
a  small  systematic  error  in  the  data.  Together  with  the  high 
values  obtained  for  the  correlations  to  logarithmic  profiles, 
these  data  indicate  an  overall  system  accuracy  of  considerably 
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better  than  10%  in  each  case  and  better  than  5%  in  perhaps 
half  of  the  data. 

tty  next  concern  was  to  examine  the  charts  to  determine 
whether  or  not  bursts  were  present,  and  whether  they  appeared 
randomly,  sequentially,  or  simultaneously  at  the  various 
levels.  Once  I  knew  the  nature  of  the  phenomenon  I  could 
begin  computing  averages  and  profile  parameters. 

The  purity  of  the  results  was  unexpected.  Each  of  the 
three  high  frequency  records  from  the  22  acceptable  runs  was 
searched  and  marked  for  bursts.  In  no  run  were  there  any 
bursts  which  did  not  occur  simultaneously  (±1  second  at  the 
leading  and  trailing  edges)  at  all  three  levels!  This  may 
well  mean  that  there  is  an  unexpectedly  rapid  exchange  of 
turbulent  properties  in  the  flow.  This  will  be  discussed 
in  subsequent  paragraphs.  The  intermittency  factor,  I, 
defined  as  the  ratio  of  time  occupied  by  bursts  to  the  length 
of  the  record,  ranged  from  0.27  to  0.82  in  an  irregular  way. 
Table  11  shows  the  distribution  of  values.  Clearly,  most 
of  them  fall  between  0.50  and  0.69-  Thus  it  might  have  been 
more  orthodox  to  have  defined  the  regions  of  lower  level 
activity  as  "bursts"  since  they  comprise  the  smaller  part 
of  the  records . 

The  tine  between  bursts  and  burst  length  seemed  to 
exhibit  no  detectable  pattern  even  though  other  experimenters 
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have  noted  that  the  burst-  rate  in  laboratory  work  correlates 
somewhat  with  friction  velocity,  Kim,  Kline,  and  Reynolds 
(1971),  at  smaller  friction  velocities  than  occurred  in  our 
work.  Typical  bursts  lasted  from  5  to  50  seconds  with  modal 
values  around  15  seconds.  The  intervals  between  bursts  were 
similar.  We  did  no  detailed  analysis  on  the  burst  frequencies 
and  lengths  after  determining  that  no  significant  correlation 
with  other  variables  was  to  be  found. 

There  were  88  profiles  resulting  from  the  22  runs.  Mean 
velocities  were  computed  from  the  hob-film  systems  for  the 
portions  of  the  records  within  bursts,  without  bursts,  and 
both  combined.  Additionally,  this  last  was  also  computed 
with  the  cup  anemometer  data  included.  The  data  were  fitted  to 

kU/u* 

z  =  z  e 
o 

The  correlation  coefficients  were  quite  high.  Table  12  shows 
their  distribution. 

The  major  factor  contributing  to  a  poor  correlation 
coefficient  at.  the  lower  levels  in  some  runs  was  the  effect 
of  lapse  rate  on  the  profile.  In  all  but  four  of  the  33 
cases  where  the  correlation  coefficient  was  smaller  than 
O.98  the  lapse  rate  in  the  lowest  meter  was  greater  in 
magnitude  than  0.5  C  per  meter  and  negative  thus  indicating 


TO 


Range  of 
Corr  Coeff* 

Cups  &  Films 
Total  Record 

Films 

Whole  Record 

Films 

Bursts 

Films 

Bursts  Elim 

TOTAL 

O.98-I.OO 

22 

10 

11 

12 

55 

0.96-0,98 

0 

6 

2 

5 

13 

O.9U-O.96 

0 

U 

6 

1 

11 

0.92-0. 9H 

0 

1 

2 

2 

5 

0.90-0.92 

0 

1 

1 

1 

3 

0.88-0.90 

0 

0 

0 

0 

0 

0 . 86-0 . 88 

0 

0 

0 

0 

0 

0.8L-0.86 

0 

0 

0 

0 

0 

0.82-0.81* 

0 

0 

0 

1 

1 

TOTAL 

22 

22 

22 

22 

88 

*  Class  intervals  are  open  to  the  right 


Table  12.  Distribution  of  correlation  coefficients 


that  strongly  unstable  conditions  obtained.  On  the  other 
hand,  only  four  correlations  less  than  0.96  occurred  with 
lapse  rates  closer  to  neutral  than  -0.5  C/in. 

The  lack  of  large  variation  between  the  goodness  of 
fit  to  the  logarithmic  formula  among  the  three  kinds  of  pro¬ 
files  made  with  the  films  may  indicate  that  the  adjustment  of 
the  profile  to  new  conditions  is  quite  rapid.  This  is  further 
supported  by  some  tests  conducted  low  over  the  furrows.  The 

test  was  intended  to  examine  the  structure  of  the  flow  over 

,  \ 

the  furrows  by  measuring  profiles  as  a  function  of  phase  with 

respect  to  the  furrows.  This  v(as  not  done  because  the  cup 
anemometers  were  not  available  on  the  day  allotted  to  the 
experiment,  but  hot-film  data  were  obtained  along  a  line  nor¬ 
mal  to  the  crest  lines  10  cm  above  the  top  of  the  crests,  see 

Fig.  23.  The  results  showed  that  the  velocity  at  10  cm  ( 1/9 

i 

wavelength,  1  peak  to  peak  amplitude)  above  the  furrows  was, 
within  ±15$,  the  same  as  if  a  logarithmic  profile  having  the 
same  roughness  length  as  obtained  for  our  other  records  had 
been  mechanically  lifted  without  change  in  shape  t<.  follow 
the  contour  of  the  furrow.  For  a  1  m/s  wind  at  probe  height 
the  time  to  pass  the  22  cm  between  probes  is  about  200  ms. 

The  data  are  inadequate  to  determine  whether  the  quick  ad¬ 
justment  to  the  change  in  position  of  the  surface  attenuates 
with  height.  The  conventional  wisdom  is  that  rapid  attenu¬ 
ation  should  occur  but  the  simultaneity  of  the  appearance  of 
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bursts  implies  so  high  a  rate  of  vertical  transport  of  flow 
properties  that  the  attenuation  may  not  occur.  The  adjustment 
of  the  profile  may  explain  why  Takeda  ( 1963 )  had  difficulty 
following  the  evolution  of  wave  induced  kinks  in  profiles 
passing  over  sand  bars,  the  kinks  disappearing  so  swiftly 
that  they  were  undetectable  by  the  time  they  reached  his 
land  stations. 

Since  the  bursts  occurred  simultaneously  at  all  levels 
and  since  the  profiles  adjusted  rapidly  enough  that  good  cor¬ 
relations  were  obtained  within  and  without  the  bursts,  our 
data  were  appropriate  for  examining  the  parameters  of  the 
logarithmic  profiles  to  test  Ruggles'  hypothesis. 

Before  I  began  the  search  for  systematic  variations 

of  z  and  u*  as  a  function  of  U  and  the  bursts ,  I  marked 
°  10 

the  22  runs  with  "flags"  indicating  in  each  of  the  4  profiles 
of  each  run  such  anomalies  as  one  would  not  expect  to  find 
in  "perfect"  data.  Flags  were  given  for  the  following  reasons 

Lapse  rate  steeper  than  -0.5  C/m 
Correlation 

u*  (burst) _ 

u*  (3  level  combined) 

u*  (burst  absent) _ 

u*  (3  level  combined) 

zQ  (6  level  combined) 
z_  (3  level  combined) 

O 


=0.98 


<1.0 


>1.0 


>2.0  or  <0.5 


Ik 


U10  (6  level  combined) 
Uio  (3  level  combined) 


>1.1  or  <0.9 


(burst  absent) 

U1Q  (3  level  combined) 


>1.0 


U^q  (burst) 

U1Q  (3  level  combined) 


<1.0 


The  largest  number  of  flags  which  could  be  placed  on  any  single 
non  was  eleven,  since  a  correlation  flag  was  given  if  one  pro¬ 
file  was  insufficiently  correlated,  two  if  two,  et  cetera. 

The  runs  were  then  classified  in  order  of  reliability  by  the 
number  of  flags  they  were  assigned,  Te.ble  13.  The  data  were 
analyzed  using  the  best  9,  best  12,  best  IT,  and  the  entire 
set  of  22. 

Because  we  wanted  to  plot  z„  as  a  function  of  U, 

e  o  10 

the  mean  wind  at  10  meters,  as  Ruggles  did,  and  because  we 
had  no  anemometer  so  high,  we  computed  U1  ^  from  the  profile 
parameters  of  the  profile  for  which  U  was  desired.  The 
result  for  all  88  profiles  is  shown  in  Fig,  2h.  Except  for 
a  hint  of  a  slight  increase  in  zQ  near  =  6.5  m/s  which 
may,  if  real,  be  due  to  some  sort  of  instability,  the  rough¬ 
ness  length,  z  ,  seemed  relatively  insensitive  to  everything, 
o 

For  each  run  I  computed  the  ratios 


Number  of  Flags  N 


<2 


3 


Number  of  Runs 
having  N  Flags 

Number  of  Runs  having 
N  or  Fewer  Flags 


9  3 

9  12 


Table  13.  Distribution  of  flags 
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zq  (burst) 

zQ  (3  level  combined) 


and 


zQ  (burst  absent) 
z0  (3  level  combined) 


and  calculated  their  geometric  and  arithmetic  means  and  standard 
deviations,  Table  lU.  The  mean  values  of  the  ratios  are  close 
to  unity  within  the  precision  of  the  system.  The  arithmetic 
and  geometric  means  and  standard  deviations  do  not  differ 
greatly.  Thus  a  separate  computation  of  the  geometric  mean 
and  standard  deviation  for  zQ  seems  unnecessary.  The  arith¬ 
metic  mean  values  and  standard  deviations  for  zQ  are  presented 
in  Table  15.  The  only  trend  i.s  a  slight  decrease  in  zQ  as  the 
quality  of  the  data  decreases.  This  happens  to  correlate  with 
an  increasingly  unstable  lapse  rate  but  the  magnitude  of  the 
trend  and  the  sample  size  are  both  so  small  as  to  be  insuf¬ 
ficient  for  drawing  any  conclusions. 

If  the  roughness  length  seems  relatively  unaffected  by 
changes  in  both  the  mean  wind  and  the  burst  structure,  what 
about  the  friction  velocity?  Figure  25  shows  u*  va.  II  for 
all  88  profiles.  Except  for  a  hint  of  scatter  near  6.5  ro/s , 
the  data  clump  uniformly  about  a  straight  line  u*  =  -0.035  + 
O.Oii69  U  ^  with  a  correlation  of  O.895.  Table  16  shows  the 
variation  of  u*  with  the  burst  structure.  While  the  magni¬ 
tude  of  the  ratios  is  of  the  same  order  as  that  for  the  ratios 

of  z  's,  the  value  for  a  is  smaller.  Thus,  these  variations 
0 
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Ail  ; 

1.100 
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0.39*4 

0.1*18 

0.385 

1.021 

0.858 

0.896 

0.406 

0.581 

0.61*1 

1.500 
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1.897 

1.136 

1.170 

1.389 

0.522 

0.600 

0.911 

O.98J 
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0.605 

0.65** 

0.681 

1.831 

1.923 

1.977 
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12  Best 


0.00127 

0.0014 

0.00123 

0.0011 

0.00104 

0.00083 


17  Best 


0.00115 

0.0012 

0.00133 

0.0012 

0.00105 

0.00075 


All  22 


0.000989 

0.0011 

0.00129 

0.0011 

0.000937 

0.00071 


devi ations 


Selected  Data 


9  Best  12  Best  IT  Best  All  22 


Bursts 

u*/u*  (3  level  combined) 

CT 

l.lUli 

0.073 

1.134 

0.073 

1.137 

0.068 

1.136 

0.108 

Bursts  Absent 

u*/u*  (3  level  combined) 

0.835 

0.830 

0.865 

0.885 

a 

0.078 

0.080 

0.092 

0.100 

Table  1 6.  u*  ratio  means  and  standard  deviations 
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are  more  significant  numerically  and  much  more  significant 
physically  since  they  enter  the  profile  linearly  while  the 
roughness  length,  zQ,  enters  the  profile  logarithmically . 

For  each  hot-film  profile  I  also  computed  the  ratio  in 
constant  arbitrary  units  of  u^/E^q  where  E^  is  a  quantity 
associated  with  the  mean  turbulent  energy  above  20  Hz.  The 
averages  of  these  values  are  presented  in  Table  IT.  The 
results  indicate  that  the  energy  in  the  filtered  portion  of 
the  turbulence  above  20  Hz  is  an  appropriate  indicator  but 
not  an  exact  measure  of  the  value  of  the  friction  velocity. 

If  ve  use  the  friction  velocity  as  the  mere  reliable  indi¬ 
cator  of  total  turbulence  activity  in  the  flow,  then  from 
Table  16  we  may  compute  che  ratio 


u*  burst _ 

u*  burst  absent 

for  the  better  data.  Thus  we  find  about  twice  as  much  turbu¬ 
lent-  energy  in  the  hursts  as  outside  of  them.  We  are  not 
sure  why  this  seems  to  occur  with  such  uniformity.  Because 
the  disparity  in  Reynolds  numbers  is  so  great  between  lab¬ 
oratory  experiments  and  those  in  the  atmosphere,  laboratory 
experiments  in  intermittent  flows  may  provide  little  information 
of  use  in  explaining  our  observations. 


83 


Selected  Data 
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a 
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3  Level  Combined 

u*2/& 

20 

0.00131 

0.000861 

a 

0.00109 

0.000940 

20 


Tabic  17.  Means  and  standard  deviations  of  u*2/E. 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  occurrence  of  bursts  simultaneously  at  all  levels  and  the 
validity  of  the  logarithmic  profile  within  and  without  them 
indicate  large  vertical  transport  of  mean  momentum  through  the 
agency  of  nearly  instantaneous  adjustment  of  the  Reynolds 
stresses.  If  the  bursts  were  significantly  swot  by  or  frozen 
into  the  mean  flow,  the  upper  levels  shoulvi  see  t  ••  hursts 
first.  That  they  do  not  is  a  good  indication  that  the  bursts 
don't  attach  to  specific  fluid  elements — they  are  purely 
dynamic  occurrences .  If  this  is  the  case,  Taylor's  hypoth¬ 
esis  for  converting  temporal  spectra  to  spatial  spectra  may 
be  inapplicable  in  the  lower  few  meters  of  the  atmospheric 
boundary  layer. 

The  lack  of  large  peaks  in  the  zQ  values  for  certain 
values  of  the  wind  at  anemometer  height,  and  the  associated 
smoothness  of  the  u*  vs.  curve  indicate  that  the  mechanism 

causing  the  peaks  in  Ruggles'  data  was  not  active  in  our 
experiment.  Since,  except  for  the  rigidity  of  the  surface, 
we  duplicated  Ruggles'  work  closely,  it  seems  reasonable  to 
conclude  that  the  mechanism  was  due  to  the  flexing  of  the 
surface.  That,  of  course,  means  trie  mechanism  was  due  to 
the  presence  of  waves  on  the  surface  and  thereby  it  is  indi¬ 
cated  that  the  waves  influence  not  merely  pen iodic  perturbations 


65 


on  the  profile,  but  the  profile  proper. 

Since  the  profile  is  a  function  of  the  state  of  the  sea, 
if  waves  are  being  generated  then  the  profile  will  be.  changing 
e.s  well  and  must  be  obtained  from  an  adequate  model  as  part  of 
the  solution.  This  complicates  the  already  difficult  problem 
of  devising  appropriate  averaged  field  equations.  The  close 
coupling  between  sea  and  air  at  the  interface  may,  however, 
simplify  the  establishment  of  boundary  conditions. 

Because  the  wind  field  contains  sharply  defined  regions 
of  higher  turbulence  activity  and  higher  horizontal  mean  velo¬ 
city  than  surrounding  regions,  the  flow  certainly  can't  be 
horizontal  and  two  dimensional.  The  equation  of  continuity 
requires  either  significant  vertical  flow  components  or  sig¬ 
nificant  variability  in  the  direction  of  the  horizontal  wind, 
or  both,  throughout  the  flow.  Many  of  our  simplifying  assump¬ 
tions  for  the  mathematical  modeling  of  the  atmospheric  boundary 
layer  axe  invalid.  Only  a  full  three  dimensional  treatment  is 
adequate  for  the  micro-climatological  scale  analysis  of  geo¬ 
physical  flows. 

The  three  dimensional  structure  may  account  for  the 
arcade  patterns  noted  on  the  sea  surface,  Kinsman  (1965), 
page  5^3.  The  turbulent  energy  within  the  bursts  is  nearly 
twice  that  outside  and  if  the  remaining  Reynolds  stresses 
scale  proportionately,  the  surface  stress  would  be  roughly 
doubled  in  the  high  activity  regions.  The  result  would  be 
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alternating  regions  of  high  and  low  wave  generation  activity 
behaving  somewhat  as  microscale  "storms"  each  radiating  its 
own  swell. 

To  make  further  progress  in  understanding  the  actual 
structure  of  the  lower  10  meters  of  the  atmosphere  over  the 
ocean  and  over  land  an  intense  investigation  of  the  three 
dimensional  structure  of  the  flow  should  be  conducted.  An 
ideal  program  would  involve  simultaneous  measurements  sensing 
3-component  velocity  profiles  over  heights  from  3  centimeters 
to  10  meters  at  each  station  and  taking  careful  records  of 
temperature  profiles  as  well  along  a  line  leading  down-wind 
from  land,  across  a  beach,  and  out  to  sea  with  the  developing 
wave  field.  The  instruments  should  be  of  fast  response  type 
to  allow  Reynolds  stresses  and  heat  fluxes  to  be  computed  to 
high  accuracy. 

This  would  be  expensive,  but  not  nearly  so  expensive 
as  continuing  the  fiction  of  a  horizontal  two  dimensional  flow. 
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APPENDIX 


Data  Compiled  from  +.he  Runs  Used  i 


Run:  1  Date:  20  Mar  71  Tape:  2  Section:  1 

Run  Length:  16.5  min 

Temperatures:  li . 5  C  at  1.5  m  and  4.85  C  at  0.5  m 
Number  of  Flags:  1  Intermittency :  0.6l 

Six  Level  Combined  Data 


z  m 

4  2 

1 

0.5  0.25 

0.125 

U  m/s 

8.76  7.98 

7.22 

6.39  5-93 

4.92 

u»:  0. 

434  m/s  zQ:  0. 

00126  m  Corr.  coeff . :  0.997  U  :  9*739  m/s 

Three 

Level  Combined 

Data 

z  m 

0.5  0.25 

0.125 

zc:  0.00141  m  u*; 

0.444  m/s 

U  m/s 

6.39  5-93 

4.92 

Corr.  coeff.:  0.977 

U1Q: 

9-  84l  m/s 

u*2/E20:  0.0011 

Three 

Level  With  Bursts  Only 

z  m 

0.5  0.25 

0.125 

zQ:  0.00176  m  u*: 

0.522  m/s 

U  m/s 

7.26  6.63 

5.50 

Corr.  coeff.:  0.987 

C 

H 

O 

11.271  m/a 

u*2/E2Q:  0.0011 

Three 

Level  With  Bursts  Omitted 

z  m 

0.5  0.25 

0.125 

zQ:  0.00202  m  u*: 

0.394  m/s 

U  m/s 

5-38  4.83 

4.03 

Corr.  coeff.:  0.994 

U  : 

10 

8. 382  m/s 

u*2/E20:  0.0029 
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Run:  2 


Date:  20  Mar  71 


Tape:  2  Section:  2 


Run  Length:  16 . ?  min 

Temperatures:  3. 85  C  at  1.5  m  and  3.9  C  at  0.5  m 
Number  of  Flags:  none  Intennittency :  0.70 


Six  Level  Combined  Data 


Z  ffi 

)l 

2 

1 

0.5  0.25 

0.125 

U  m/s 

7  •  54 

6.84 

6,22 

5.49  4.91 

4.23 

u*:  0. 

381  m/s  zQ:  0.00148  m  Corr.  coei'f . :  1.000  U^q:  8.39  m/s 

Three 

Level 

Combined 

Data 

z  m 

0.5 

0.25 

0.125 

zQ:  0.00118  m  u* : 

0.364  m/s 

U  m/s 

5.49 

4.91 

4.23 

Corr.  coeff . :  1.000 

U10  = 

8.238 

m/s 

u*2/E2Q:  0.0014 

Three 

Level 

With  Bursts  Only 

z  m 

0.5 

0.25 

0.125 

zQ:  0.00164  m  u*: 

0.425  m/s 

U  m/s 

6. 06 

5.38 

4.59 

Corr.  coeff:  0.999 

U10: 

9.262 

m/s 

u*2/E2Q:  0.0015 

Three 

Level 

With  Bursts  Omitted 

z  m 

0.5 

0.25 

0.125 

zQ:  0.00237  m  u#: 

0.336  m/s 

U  m/s 

4.53 

3.84 

3.38 

Corr.  coeff. :  0.993 

U10: 

7-015 

m/s 

u*2/E2q:  0.0036 
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Run:  3  Date:  25  Mar  71 
Run  Length:  10.5  min 
Temperatures:  0.75  0  at  1.5  m 
Number  of  Flags :  none 


Tape:  1  Section:  1 

and  1.25  C  at  0,5  m 
Intermittency :  0.63 


Six  Level  Combined  Data 


z  m 

4 

2  1 

U  m/s 

5.53 

4.96  4.47 

u*:  0. 

290  m / 

's  z  :  0.00200  m  Cc 
0 

Three 

Level 

Combined  Data 

z  m 

0.5 

0.25  0.125 

U  m/s 

4.01 

3.64  2.92 

U10: 

6.522 

m/s 

Three 

Level 

With  Bursts  Only 

z  m 

0.5 

0.25  0.125 

U  m/s 

4.60 

3.98  3.20 

U10: 

7.669 

m/s 

Three 

Level 

With  Bursts  Omitted 

z  m 

0.5 

0.25  0.125 

U  m/3 

3.40 

3.05  2.43 

uio: 

5.608 

m/  s 

0.5  0.25  0.125 

4.01  3.64  2.92 

Corr.  coeff.:  0.997  U^:  6.176  m/s 

zq:  0.00328  m  u#:  0.325  m/s 

Corr.  coeff.:  O.983 
u*2/E20:  0.0038 

0.00520  m  u*:  0.406  m/s 
o 

Corr.  coeff.:  0.99® 
u*2/E20:  0.0046 

z0:  0.00405  m  u*:  0.287  m/s 
Corr.  coeff.:  O.987 
u*2/E20:  0.0050 
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Run:  4 


Date:  25  Mar  71  Tape:  1  Section:  2 


Run  Length:  10.75  min 

Temperatures:  0.9  C  at  1.5  m  and  1.6  C  at  0.5  m 
Number  of  Flags:  1  Intermittency :  0.67 

Six  Level  Combined  Data 


z  m 

4 

2  1 

0.5  0.25 

0.125 

U  m/s 

5.60 

5.04  4.55 

4.09  3.68 

3-13 

u*:  0. 

,279  m/a  zQ:  0.00139  m  Corr.  coeff . :  0.999  upo:  6.198  m/s 

Three 

Level 

Combined  Data 

z  m 

0.5 

0.25  0.125 

zQ:  0.00137  m  u*: 

0.279  m/s 

U  m/s 

4.09 

3.68  3.13 

Corr.  coeff.:  0.996 

uio: 

6.203 

m/s 

u*2/E20i  0.0016 

Three 

Level 

With  Bursts  Only 

z  m 

0.5 

0.25  0.125 

zq:  0.00215  m  u*: 

0.346  m/s 

U  m/s 

4.72 

4.11  3.52 

Corr.  coeff.:  1.000 

U10: 

7-312 

m/s 

u*z/E„„:  0.0018 

20 

Three 

Level 

With  Bursts  Omitted 

z  m 

0.5 

0.25  0.125 

zQ:  0.00148  m  u*: 

0.214  m/s 

U  m/s 

3.08 

2.79  2.35 

Corr.  coeff.:  0.993 

U10: 

4.708 

m/s 

u»2/E20:  0.0046 
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Run: 


Date:  25  Mar  71  Tape:  1  Section:  3 


Run  Length:  10.25  min 

Temperatures:  1.5  C  at  1.5  m  and  2.00  C  at  0.5  m 
Number  of  Flags:  1  Intermit tency :  0.60 

Six  Level  Combined  Data 


z  m 

It 

2 

1 

0-5 

0.25 

0.125 

U  m/s 

5.18 

4.65 

4.16 

3.70 

3.40 

2.88 

u*:  0. 

260  m/s 

zQ:  0.00153  m 

Corr.  coeff.: 

0.997 

U1Q:  5.719  m/s 

Three  Level  Combined  Data 


z  m 

0.5 

0.25 

0.125 

zQ:  0.00103  m  u" :  0.242  m/s 

U  m/s 

3.70 

3.40 

2.88 

Corr.  coeff.:  0.988 

U10: 

5.561 

m/s 

u*2/E20:  0.0033 

Three 

Level 

With  Bursts 

Only 

z  m 

0.5 

0.25 

0.125 

zQ :  0.000864  m  u*:  0.265  m/s 

U  m/s 

4.15 

3.84 

3.26 

Corr.  coeff.:  0.985 

U10: 

6.192 

m/s 

vi*’/E20:  0.0028 

Three 

Level 

With  Bursts 

1  Omitted 

z  m 

0.5 

0.25 

0.125 

zc:  0.000747  m  u*:  0.l8l  m/s 

U  m/s 

2.86 

2.75 

2.30 

Corr.  coeff.:  0.944 

O 

1 — 1 

4.309 

m/s 

u»2/E20:  0.0044 

lib 


\ 
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Run:  6  Date.:  8  April  ;71 

Run  Length:  10.0  min 

Temperatures:  7*8  C  at  1.0  m 

Number  of  Flags :  5 

\ 

i 

Six  Level  Combined  Data 


z  m 

4 

2 

1 

;U  m/s 

7.66 

6.81 

6.07 

u*:  0. 

377  m/s  zQ:  0 

.00135  m 

Three 

Level 

Combined 

Data  \ 

z  m 

0.25 

0.12 

0.06 

U  m/s 

4.78 

4.48 

3.55 

uio: 

8.418 

\ 

m/s 

\ 

Three 

Level 

With  Bursts  Only 

z  m 

0.25 

0.12 

0.06 

LI  m/s 

4.94 

4.68 

3.75 

U10: 

8.543 

m/s 

Three 

Level 

With  Bursts  Omitted 

z  m 

0.25 

0.12 

o'.  06 

U  m/s 

3.48 

3.53 

\ 

2.66 

U10: 

6.854 

m/s 

Tape:  1  Section:  1 

\ 

\ 

1  I 

and  8.5  C  at  0.2  m 
Intermittency :  0.82 

0,25  0.12  0.06 

4.78  4.18  3-55 

» 

i 

Corr.  coeff.!  0.995  U10:  8.390  m/s 

zQ:  0.00131  m  u*:  0.377  m/s 

Corr.  coeff.:  0.954 

U*2/E20:  °- 000141 

l 

1 

z0:  0.000993  m  u*:  0.371  m/s 

Corr.  coeff.:  0.946 
u*2/E2Q:  0.00033 

zQ:  0.00245  m  u*:  0.330  m/s 

Corr.  coeff.:  0.83 
u*2/E20:  0.0019-. 
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Run:  7  Date:  8  April  71  Tape:  1  Section:  2 

Run  Length:  11.0  min 

Temperatures:  7*9  C  at  1.0  m  and  8.6  C  at  0.2  m 
Number  of  Flags:  4  Intermittency:  0.39 

Six  Level  Combined  Data 


z  m 

4 

2 

1 

0.25  0.12 

0.06 

U'  m/s 

7.07 

6.35 

5.68 

4.50  4.15 

3.28 

u*:  0. 348  m/s 

v  °- 

1 00129  m  Corr.  coeff.:  0.997  U10:  7-788  m/s 

\ 

Three 

Level  Combined 

i 

Data 

\ 

z  m 

0.25 

0.12 

0.06 

i 

zQ:  0.00154  m  u*: 

0.364  m/s 

U  m/s 

4.50 

4.15 

3.28 

Corr.  coeff . :  O.967 

U10: 

7.995  m/s 

u*2/E20:  0.00060 

Three 

Level  With  Bursts  Only 

\ 

\ 

\ 

0.427  m/s 

z  m 

0.25 

0.12 

0.06 

zq:  0.00155  m  u*: 

U  m/s 

5.21 

4.91 

3.84 

Corr.  coeff.:  0.946 

U10! 

9-362  m/ 

s 

u*2/E20:  0.00044 

Three 

Level  With  Bursts  Omitted 

z  m 

0.25 

0.12 

0.06 

zQ:  0.00143  m  u*: 

0.318  m/s 

U  m/s 

4.02 

3.66 

2.93 

Corr.  coeff.:  0.978 

V 

7.043  m/ 

s 

u#2/E2():  0.00099 
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Run:  8  Date:  8  April  71  Tape:  1  Section:  3 

Run  Length:  11.0  min 

Temperatures:  8.1  C  at  1.0  m  and  8.85  C  at  0.2  m 
Number  of  Flags:  1+  Intemittency :  O.56 

Six  Level  Combined  Data 

z  m  1*  2  1  0.25  0.12  0.06 

U  m/s  8.79  7.81  6. 94  5.69  5-52  4.36 

u*:  0.399  m/s  zQ:  0.000716  m  Corr.  coeff.:  O.988  U1Q:  9-521  m/s 

Three  Level  Combined  Data 

z  m  0.25  0.12  0.06  zQ'.  0.00114  m  u*:  0.445  m/s 

U  m/s  5.69  5.52  4.36  Corr.  coeff.:  0.912 

U1Q:  10.099  m/s  u*»2/E20:  0.00022 

Three  Level  With  Bursts  Only 

z  m  0.25  0.12  0.06  zQ:  0.00112  m  u*:  0.1+91  m/s 

U  m/s  6J+0  6.05  4.82  Corr.  coeff.:  0.947 

U1Q:  11.176  m/s  u*2/E20:  0.00019 

Three  Level  With  Bursts  Omitted 

z  m  0.25  0.12  0.0 6  z0:  0.00149  m  u*:  0.4l4  m/s 

U  m/s  5.08  4.82  3-77  Corr.  coeff.:  0.939 

U10:  9.121  m/s  u*2/E2C:  0.00040 


99 


Run:  9  Date:  8  April  71 

Run  Length:  10.5  min 
Temperatures:  9*2  C  at  1.0  m 
Number  of  Flags:  3 

Six  Level  Combined  Data 
z  m  4  2  1 

U  m/s  8.05  7-25  6.50 

u*:  0.328  m/s  zq:  0.000287  m 

Three  Level  Combined  Data 
z  m  0.25  0.06  0.03 

U  m/s  5.14  4.49  4.10 

Uj^q :  6.948  m/s 

Three  Level  With  Bursts  Only 
z  m  0.25  0.06  0.03 

U  m/s  5-69  4.99  4.64 

U10:  7.518  m/s 

Three  Level  With  Bursts  Omitted 
z  m  0.25  0.0b  0.03 

U  m/s  4.07  3.71  3.35 

U1Q:  5.359  m/s 


Tape:  2  Section:  1 

and  9.85  C  at  0.2m 
Intermittency :  0.59 

0.25  0.06  0.03 

5.14  4.49  4.10 

Corr.  coeff . :  0.986  U-^q!  8.58  m/s 

zc-  0.0000063  m  u*:  0.195  m/s 

Corr.  coeff.:  0.998 
u*2/E20:  0.00007 

zQ:  0.0000025  m  u*:  0.198  m/s 

Corr.  coeff.:  1.000 
u*2/E20:  0.000048 

zQ:  0.0000014  m  u*:  0.136  m/s 

Corr.  coeff.:  0.98l 
u*2/E20:  0.00010 


»fiu 


Run:  10  Date:  8  April  71  Tape:  2  Section:  2 
Run  Length:  11.0  min 

Temperatures:  10.2  C  at  1.0  m  and  10.8  C  at  0.2  m 
Number  of  Flags:  4  Intermittency :  0.50 

Six  Level  Combined  Data 


z  m 

1+  2 

1 

0.25  0.06 

0.03 

U  m/s 

6.74  6.  04 

5.39 

4.27  3.49 

3.24 

u*:  0.292  m/s  z0:  0.000499  m  Corr.  coeff . :  0.991  U^q:  7.225  m/B 

Three 

Level  Combined 

Data 

z  m 

0.25  0.06 

0.03 

zQ:  0.000494  m  u*: 

0.200  m/s 

U  m/s 

4.27  3.^9 

3.24 

Corr.  coeff.:  0.9 96 

U10: 

4.948  m/s 

u*2/E2Q:  0.00020 

Three 

Level  With  Bursts  Only 

z  m 

0.25  0.06 

0.03 

z0:  0.000723  m  u*: 

0.238  m/s 

U  m/s 

4.85  3-97 

3.60 

Corr.  coeff.:  0.999 

UlO* 

5.665  m/s 

u*2/E20:  0.00018 

Three 

Level  With  Bursts  Omitted 

z  m 

0.25  0.06 

0.03 

z0:  0.000126  m  u*: 

0.203  m/s 

U  m/3 

3.87  3-02 

2.87 

Corr.  coeff. :  O.983 

U10: 

5.727  m/s 

u*2/E20:  0.00044 

Run:  11  Date:  8  April  71  Tape:  2  Section:  3 
Run  Length:  10.5  min 

Temperatures:  10.5  C  at  1.0  m  and  11  C  at  0.2  m 
Number  of  Flags:  3  Intermittency :  0.77 

Six  Level  Combined  Data 


z  m 

4  2 

1 

0.25  0.06 

0.03 

U  m/s 

7.88  7-04 

6.31 

4.92  4.07 

3.81 

u*:  0.342  m/s  zQ:  0. 

,000515  m  Corr.  coeff.:  0.988 

U1Q:  8.439  m/s 

Three 

Level  Combined 

Data 

z  m 

0.25  0.06 

0.03 

zQ:  0.0000283  m 

u* :  0 . 21 6  m/  s 

U  m/s 

4.92  4.07 

3.81 

Corr.  coeff.:  0. 

995 

U10: 

6.894  m/s 

u*2/E2Q:  0.0010 

Three 

Level  With  Bursts  Only 

z  m 

0.25  0.06 

0.03 

zQ:  0.0000487  m 

u*:  0.247  m/'s 

U  m/s 

5.30  4.31 

4.04 

Corr.  coeff. :  0. 

993 

u10: 

7.560  m/s 

u*2/®20:  0-00011 

Three 

Level  With  Bursts  Omitted 

z  m 

0.25  0.06 

0.03 

zQ:  0.0000559  m 

u*:  0.190  m/s 

U  m/s 

4.00  3.26 

3.02 

Corr.  coeff,:  0. 

996 

U10: 

5-737  m/s 

U*2/E20:  °-00025 

Run:  12 


Date:  8  April  71  Tape:  3  Section:  1 


Run  Length:  10.25  min 

Temperatures:  10.1:  C  at  1.0  m  and  11.0  C  at  0.2  a 
Number  of  Flags:  2  Intermittency :  0.27 


Six  Level  Combined  Data 


z  m 

l * 

2 

1 

0.25  0.06 

0. 

U  m/s 

7.81 

7.01+ 

6.30 

5.02  3.95 

3. 

u*:  0 

.  367  m/ s 

z0:  0. 

000909  m 

Corr.  coeff, :  0.998  U 

10 : 

Three 

Level  Combined 

Data 

7,  m 

0.25 

0.06 

0.03 

zQ:  0.000592  m  u*: 

0 

U  m/s 

5.02 

3-95 

3.23 

Corr.  coeff.:  0.996 

u10: 

8.139  m/s 

u*2/E20:  0.00032 

Three 

Level  With  Bursts  Only 

z  m 

0.25 

0.06 

0.03 

zQ:  0.000452  m  u*: 

0 

U  m/s 

5-94 

4.96 

3.91 

Corr.  coeff.:  0.977 

u10! 

9.620  m/s 

u*2/E20:  0.00018 

Three 

Level  With  Bursts  Omitted 

z  m 

0.25 

0.06 

0.03 

zQ;  0.000737  m  u*: 

0 

U  m/s 

4.69 

3.55 

2.98 

Corr.  coeff.:  1.000 

U10: 

7.66l  m/s 

u*2/E20:  0.00075 

03 

23 

8.530  m/s 

.  331*  m/s 

.385  m/s 

.322  m/s 


Run:  13 


Date:  8  April  71 
Run  Length:  10.5  min 
Temperatures :  11.0  C  at  1.0  m 

Number  of  Flags :  1 

Six  Level  Combined  Data 


z  m 

4 

2 

1 

U  m/s 

7.02 

6.28 

5.58 

u*:  0.354  m/s  zQ: 

0.00161  m 

Three 

Level 

Combined 

Data 

z  m 

0.25 

0.06 

0.03 

U  a/s 

4.45 

3.24 

2.61 

U10! 

7.646 

m/s 

Three 

Level 

With  Bursts  Only 

z  m 

0.25 

0.0  6 

0.03 

U  m/s 

4.97 

3.67 

2.95 

c 

H 

O 

8.481 

m/s 

Three 

Level 

With  Bursts  Omitted 

z  m 

0.25 

0.06 

0.03 

U  m/s 

3.64 

2.72 

2.21 

V 

6.126 

m/a 

Tape:  3  Section:  2 

and  11.5  C  at  0.2m 
Intermittency :  0.54 

0.25  0.06  0.03 

4.45  3.24  2.61 

Corr.  coeff. :  0.999  U1Q:  7.729  m/s 

Z0:  0.00145  m  u#:  0.346  m/s 

Corr.  coeff. :  0.999 

u*2/E20:  0.00078 

zQ:  0.00130  m  u*:  0.379  m/s 

Corr.  coeff.:  0.999 
u*2/E2Q;  0.00061 

z0:  0.00108  m  u*:  0.268  m/s 

Corr.  coeff.:  0.999 
u*2/E20:  0.0013 


'  04 


Run:  Ik 


Date:  8  April  71  Tape:  3  Section:  3 


Run  Length:  10.5  min 

Temperatures:  10.8  0  at  1.0  m  and  11.4  C  at  0.2  m 
Number  of  Flags:  3  Intermittency:  0.48 


Six  Lev  •’  •  'mbined  Date 


z  m 

1} 

2 

1 

0.25  0.06  0, 

U  m/s 

6.76 

6.05 

5.44 

4.35  3.62  2 

u*:  0.309  m/s  zQ:  0. 

,000751  m  Corr.  coeff . :  0.995  U^q 

Three 

Level 

Combined 

Data 

2  m 

0.25 

0.06 

0.03 

zQ:  0.000587  m  u*:  < 

U  m/s 

4.35 

3.62 

2.80 

Corr.  coeff.:  0.973 

U10: 

7.180 

m/s 

u*2/E20:  0.00050 

Three 

Level 

With  Bursts  Only 

z  m 

0.25 

0.06 

0.03 

z0:  0.000479  m  u*r 

U  m/s 

5.11 

4.49 

3.40 

Corr.  coeff. :  0.938 

U10: 

8.494 

m/s 

u*2/E20:  0.00039 

Three 

Level 

With  Bursts  Omitted 

z  m 

0.25 

0.06 

0.03 

zQ:  0.000815  m  u*: 

U  m/s 

3.60 

2.82 

2.24 

Corr.  coeff. :  0.994 

U10: 

5.982 

m/s 

u*2/£20:  °-0013 

03 

80 

7.339  m/s 

.295  m/s 

1.342  m/s 

1.254  m/s 


Run:  15 


Date:  8  April  71  Tape:  4  Section:  1 


Run  Length:  11.5  min 

Temperatures:  10.75  C  at  1.0  a  and  11.3  C  at  0.2  m 
Number  of  Flags:  4  Intermitt ency:  0.29 

Six  Level  Combined  Data 

z  m  4  2  1  0.25  0.06  0.03 

U  m/s  5.75  5.22  4.76  3.83  3.36  2.52 

u*:  0.250  m/s  zQ:  0.000446  m  Corr.  coeff . :  0-992  U^q:  6.265  m/s 


Three  Level  Combined  Data 
z  m  0.25  0.06  0.03 

U  m/s  3-83  3.36  2.52 

U10:  6.432  m/s 

Three  Level  With  Bursts  Only 
z  m  0.25  0.06  0.03 

U  m/s  5-09  4.43  3.48 

U-^q:  8.15].  m/s 

Three  Level  With  Bursts  Omitted 

z  m  0.25  0.06  0.03 

U  m/s  3.25  2.91  2.12 

U  :  5.626  m/s 

10 


zQ:  0.000550  m  u*:  0.262  m/s 
Corr.  coeff.:  0.936 
u*2/E2q:  0.00066 

zQ:  0.000304  m  u*:  0.313  m/s 

Corr.  coeff.:  0.955 
u*2/E2q:  0.00030 

z0:  0.000702  m  u*:  0.235  m/s 

Corr.  coeff.:  0.912 

U*2/E20:  0,0019 


» Ub 


Run:  16  Date:  8  April  71  Tape:  4  Section:  2 
Run  Length:  11.0  min 

Temperatures:  11.15  C  at  1.0  m  and  11.72  C  at  0.2  m 
Number  of  Flags:  5  Intermittency :  O.58 

Six  Level  Combined  Data 


z  m  4 

2 

1 

0.25 

0.06 

0.03 

U  m/s  7.22 

6.47 

5.81 

4.60 

4.03 

3.12 

uM:  0.32  m/s 

zQ:  0.000584  m 

Corr.  coeff. 

:  0.991  U10:  7-799  m/s 

Three  Level 

Combined 

Data 

z  m  0.25 

0.06 

0.03 

z0:  0.000357  m  u»: 

0.231  m/s 

U  m/s  4.60 

4.03 

3.12 

Corr.  coeff.:  0.946 

U10:  7-469  m/s 

u*2/E20: 

0.00019 

Three  Level  With  Bursts  Only 

2  m  0.25  0.06  0.03  zQ:  0.000359  m  u*:  0.346  m/s 

U  m/s  5.52  4.70  3.71  Corr.  coeff . :  0.9 69 

U10:  8.361  m/s  u*2/E2Q:  0.00027 

Three  Level  With  Bursts  Omitted 

z  m  0.25  0.0S  0.03  2q:  0.000534  m  u#:  0.238  m/s 

U  m/s  3.47  3-10  2.31  Corr.  coeff.:  0.920 

u*2/E20:  0.00090 


U1Q:  5.863  m/s 


Run:  17  Date:  8  April  71  Tape:  4  Section:  3 
Run  Length:  10.5  min 

Temperatures:  11.3  C  at  1.0  m  and  11.65  C  at  0.2  m 
Number  of  Flags:  2  Intermittency :  0.63 

Six  Level  Combined  Data 


z  m 

4 

2 

1 

0.25  0.06 

0.03 

l)  m/s 

7-31 

6.68 

6.09 

4.84  4.22 

3.24 

u*:  0, 

.319  m/s  z  :  0, 
0 

,000456  m 

Corr.  coeff. :  0.993  UIQ:  7-984  m/s 

Three 

Level 

Combined 

Data 

z  m 

0.25 

0.06 

0.03 

zQ:  0.000421  m  u*: 

0.315  m/s 

U  m/s 

4.84 

4.22 

3.24 

Corr.  coeff.:  0.947 

U10: 

7.933 

m/  s 

u*2/E2Q:  0.00032 

Three 

Level 

With  Bursts  Only 

z  m 

0.25 

0.06 

0.03 

zQ:  0.000506  m  u*: 

0.352  m/s 

U  m/s 

5.17 

4.63 

3.46 

Corr.  coeff.:  0.pl8 

U10: 

8.707 

m/s 

u*2/E20:  0.00029 

Three 

Level 

With  Bursts  Omitted 

z  m 

0.25 

0.06 

0.03 

zQ:  0.000189  m  u*: 

0.232  m/s 

U  m/s 

4.10 

3.49 

2.87 

Corr.  coeff.:  O.98O 

uio: 

6.316 

m/s 

u*2/E2Q:  0.00050 

Run :  18 


Date;  10  April  71  Tape:  1  Section:  1 


Run  Length:  5 .5  min 

Temperatures:  8.8  C  at  1.0  m  and  9.3  C  at  0.2  m 
Number  of  Flags:  2  Intermittency:  O.76 


Six  Level  Combined  Data 

z  m 

4  2 

1 

0.25  0.125 

0.058 

U  m/s 

9.30  8.40 

7.64 

5.94  4.95 

4. 40 

u*:  0, 

,476  m/s  zQ:  0.00167  m  Corr.  coeff . :  0.998  Upo:  19*354  m/s 

Three 

Level  Combined 

Data 

z  m 

0.25  0.125 

0.058 

z0;  0.00113  m  u»: 

0,435  m/s 

U  m/s 

5*94  4.95 

4. 4o 

Corr.  coeff.:  0.981 

Ul0  = 

9.889  m/s 

u*2/E2Q:  0.00056 

Three 

Level  With  Bursts  Only 

z  m 

0.25  0.125 

0.058 

z0:  0.000868  m  u»: 

0.435  m/s 

U  m/s 

6.22  5.25 

4.67 

Corr.  coeff,:  0.985 

u10: 

10.176  m/s 

u*2/E20:  0.00046 

Three 

Level  With  Bursts  Omitted 

z  m 

0.25  0.125 

0.058 

zQ:  0.00127  m  u*: 

0. 362  m/s 

U  m/s 

4.83  3.99 

3.56 

Corr.  coeff.:  0.977 

U10: 

8.112  m/s 

u*2/e?o:  0.0010 

Run:  1? 


Date:  10  April  71  Tape:  1  Section:  2 


Run  Length:  5*0  min 

Temperatures:  8.5  C  at  1.0  m  and  9-0  C  at  0.2  m 
Number  of  Flags:  7  Intermittency :  O.65 


Six  Level  Combined  Data 

z  n 

It 

2 

1 

0.25  0.125 

0.058 

U  m/s 

9.19 

8.39 

7.56 

5.45  4.64 

4.32 

u*:  0. 

501  m/ 

's  zQ:  0.00251*  m  Corr.  coeff.:  0.993  U1Q:  10. 361  m/s 

Three 

Level 

Combined 

Data 

z  m 

0.25 

0.125 

0.058 

zD:  O.OOO369  m  u*: 

0.331  m/s 

U  m/s 

5.^5 

4.64 

4.32 

Corr.  coeff.:  0.962 

u10: 

8.1t50 

m/s 

u*2/E2Q:  0.C0026 

Three 

Level 

With  Bursts  Only 

z  m 

0.25 

0.125 

0.058 

zQ:  0.000242  m  u*: 

0.344  m/ s 

U  m/s 

6.02 

5.15 

4.88 

Corr.  coeff.;  0.948 

u10: 

9.139 

m/s 

u*2/F20:  0.00020 

Three 

Level 

With  Bursts  Omitted 

z  m 

0.25 

0.125 

0.058 

z0:  0.00172  m  u*: 

0. 361  m/s 

U  m/s 

4.55 

3.70 

3.29 

Corr.  coeff.:  0.974 

U10: 

7.823 

m/s 

u*2/E20:  0.0010 

\ 

I 


Run:  20  Date:  10  April  71  Tape:  1  Section; 
Run  Length:  8.0  min 

Temperatures ;v  8.65  C  at  1.0  m  and  9.25  C  at  0.2  n> 

\  !  \  : 

'Number  of  Flags:  6  Intermittency:  0.5 


Six  Level  Combined  Data 


a  m 

4  2 

1 

0.25  0.125 

U  m/s 

9.56  8.64 

7-72 

6.02  5.11 

\ 

u*:  0.L7I+  m/s  z0:  C. 

00138  m  Corr.  coeff r :  0.994  U 

Three 

L^vel  Combined 

Data 

'|Z  m 

i 

0.25  0-125 

0.058 

z0:  0.000356  m  u# : 

U  m/s 

6.02  5.11 

4.80 

Corr.  coeff.:  0.953 

U10: 

9.320  m/s 

u#2/E20:  0.00032 

Three 

\ 

Level  With  Burb 

ts  Only 

' 

z  m 

0.25  0.125 

0.058 

z0:  0.00041  m  u*; 

U  m/s 

6.59  5-54 

5.27 

Corr.  coeff.:  0.936 

U10: 

10.298  m/s 

u*2/E20:  0.00030 

Three 

Level  With  Bursts  Omitted 

z  m 

0.25  0.125 

0.058 

zQ:  0.000296  m  u*: 

U  m/s 

5-22  4.51 

4.14 

Corr.  coeff.;  0.978 

c 

o 

8.003  m/s 

u*2/e2o.-  0.00047 

!U 


0.058 

4.80 

•.a :  10.  544  m/s 

\ 

0.361+  m/s 

0. Uo8  m/s 

0.307  m/s 


I 


Run:  21  Date:  10  April  71  Tape:  1  Section:  4 
Run  Length:  6.0  min  \ 

Temperatures:  8.8  C  at  1.0  m  and  9*5  C  at  0.2  m 
Number  of  Flags:  5  Intermitt ency :  0.63 


Six  Level  Combined  Data 

z  in 

4 

2 

1 

0.25  0.125 

0.058 

U  m/s 

9.01 

8.18 

7.27 

5.49  4.70 

4.30 

u#:  0.1+70  m/s  zQ :  0. 

00197  m  Corr.  coeff.:  0.995  U1Q:  10.026  m/s 

! 

Three 

Level 

Combined 

D&ta 

1 

l 

z  m  1 

0.25 

0.125 

0.058 

zc:  0.00041  m  u#s 

0.339  m/s 

U  m/s 

5.49 

4.70 

4.30 

Cor(r.  coeff.:  0.976 

\  i 

”10’ 

8.568 

m/s 

, 

u*2/E20:  0.00039 

\ 

Three 

Level 

With  Bursts  Only 

z  m 

0.25 

0.125 

O.O58 

z0:  0.000254  m  u#: 

0.342  m/s 

U  m/s 

5.92 

5.20 

4.70 

Corr.  coeff. :  0.991 

U10: 

9-043 

m/s 

u#2/E20:  0.00029 

Three 

Level 

1 

With  Bursts  Omitted 

\ 

z  m 

G.25 

0.125 

0.058 

z0:  0.000724  m  u*: 

0.317  m/s 

U  m/s 

4.69 

3.91 

3.61 

Corr.  coeff.:  O.96I 

Jio: 

7  •  566 

m/s 

u*2/E„ :  0.00083 

cl  U 

t  t  0 

*md 


Run: 

22  Date:  10  April  71 

Tape:  1  Section: 

5 

Run  Length:  7.0  min 

Temperatures:  8.9  C 

at  1.0  m 

and  $.8  C  at  0.2  m 

Number  of  Flags:  4 

Intermittency:  0.65 

Six  Level  Combined  Data 

z  m 

L  2 

1 

0.25  0.125 

0.058 

U  m/s 

7-34  6.6l 

5.97 

4.69  3.80 

3.43 

u*:  0. 381  m/s  zQ:  C 

1.00189  m 

Corr.  coeff.:  0.997  U1Q:  8.173  m/s 

Three 

Level  Combined 

Data 

2  m 

0.25  0.125 

0.058 

zQ:  0.00161  m  u*: 

O.367  m/s 

U  m/s 

4.69  3.80 

3.43 

Corr.  coeff.:  0.966 

c: 

h* 

0 

8.013  m/s 

u*2/E20:  0.00094 

Three 

Level  With  Bursts  Only 

z  m 

0.25  0.125 

0. 058 

zQ:  0.00118  m  u*: 

0.379  m/s 

U  m/s 

5.13  4.16 

3.89 

Corr.  coeff. :  0.94l 

U10: 

8.573  m/s 

u*2/E20:  0.00072 

Three 

Level  With  Bursts  Omitted 

z  m 

0.25  0.125 

0.058 

zQ :  0.00410  m  u*: 

0.377  m/s 

U  m/s 

3.91  3.15 

2.54 

Corr.  coeff.:  0-996 

Vw’- 

7.358  m/s 

u*?/Eo0:  0.0031 

